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versus  the  leading  calorimeter  jet  Px  (uncorrected)  with  |Ti(jet#l)|  <  2  for 
PYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows  the  ratio  of  the 
detector  level  to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Px  (i.e. 
response  factor) 64 

4-22  Method  1  response  factors  for  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function 
of  the  leading  jet  Px,  where  PTsum  includes  charged  particles  with  px  >  0.5 
GeV/c  and  hi  <  1  and  the  ETsum  includes  all  particles  with  |r||  <  1.  Shows  the 
particle  level  prediction  (GEN)  versus  the  leading  particle  jet  Px  and  the  detector 
level  result  (CDFSIM)  versus  the  leading  calorimeter  jet  Px  (uncorrected)  with 
|ri()et#l)|  <  2  for  PYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows 
the  ratio  of  the  detector  level  to  the  particle  level,  CDFSIM/GEN,  versus  the 
leading  jet  Px(/.e.  response  factor) 65 

4-23  Leading  jet  Px  correction  used  in  method  2  for  "leading  jet"  events.  Shows  the 
difference  in  the  observed  leading  jet  Px  at  the  detector  level  (i.e.  in  the 
calorimeter)  compared  with  the  true  Px  (i.e.  corrected)  of  matched  leading 
particle  jets  using  PYTHIA  Tune  A  and  HERWIG 66 

4-24  Method  2  response  factors  from  PYTHIA  Tune  A  for  "leading  jet"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px.  Shows  the  particle  level 
prediction  (GEN)  versus  the  leading  particle  jet  Px  and  the  detector  level  result 
(CDFSIMcor)  versus  the  leading  calorimeter  jet  Px  (corrected)  with  |T|(jet#l)|  < 
2.  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIMcor/GEN,  versus  the  leading  jet  Px  (i.e.  response  factor) 67 

4-25  Compares  the  method  1  response  factors  versus  the  leading  jet  Px  (uncorrected) 
with  the  method  2  response  factors  versus  the  leading  jet  Px  (corrected)  from 
PYTHIA  Tune  A 68 


xu 


5-1  Data  at  1.96  TeV  corrected  to  the  particle  level  using  method  1  and  method  2 
compared  with  PYTHIA  Tune  A  and  HER  WIG  at  the  particle  level.  Shows  the 
density  of  charged  particles,  dNchg/dT|d(t)  (top),  the  PTsum  density  of  charged 
particles,  dPTsum/dTid(t)  (middle),  (pr  >  0.5  GeV/c  and  |r||  <  1),  and  the  ETsum 
density,  dET/dTid(|)  (bottom),  for  particles  with  |ti|  <  1  in  the  "transverse"  region 
(average  of  "transMAX"  and  "transMIN")  for  "leading  jet"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Ft 72 

6-1  Dlustration  of  correlations  in  azimuthal  angle  A(t)  relative  to  the  direction  of  the 
leading  jet  (MidPoint,  R  =  0.7)  in  the  event,  jet#l.  The  angle  Ac])  =  (])  -  (l)jet#l  is 
the  relative  azimuthal  angle  between  charged  particles  and  the  direction  of  jet#l. 
The  "toward"  region  is  defined  by  |A(j)|  <  60"  and  JTil  <  1,  while  the  "away" 
region  is  |A(|)|  >  120°  and  |ti|  <  1.  The  "transverse"  region  is  defined  by  60°  < 
|A(t)|  <  120°  and  |r||  <  1.  Each  of  the  three  regions  "toward",  "transverse",  and 
"away"  and  has  an  overall  area  in  Tj-cj)  space  of  At|A(|)  =  4^/3.  We  examine 
charged  particles  in  the  range  pj  >  0.5  GeV/c  and  |t||  <  1,  but  allow  the  leading 
jet  to  be  in  the  region  |Ti(jet#l)|  <  2 75 

6-2  Illustration  of  correlations  in  azimuthal  angle  A(|)  relative  to  the  direction  of  the 
leading  jet  (highest  Ex  jet)  in  the  event,  jet#l .  The  angle  A(])  =  (J)  -  (|)jet#l  is  the 
relative  azimuthal  angle  between  charged  particles  and  the  direction  of  jet#l. 
The  "toward"  region  is  defined  by  \A(^\  <  60°  and  |t||  <  1,  while  the  "away" 
region  is  \A^\  >  120°  and  |ti|  <  1.  The  two  "transverse"  regions  60°  <  Acjx  120° 
and  60°  <  -Acj)  <  120°  are  referred  to  as  "transverse  1"  and  "transverse  2".  Each 
of  the  two  "transverse"  regions  have  an  area  in  r\-^  space  of  AriAcj)  =  47c/6.  The 
overall  "transverse"  region  defined  in  Fig.  3  corresponds  to  combining  the 
"transverse  1"  and  "transverse  2"  regions.  Events  in  which  there  are  no 
restrictions  placed  on  the  on  the  second  highest  Et  jet,  jet#2,  are  referred  to  as 
"leading  jet"  events  (left).  Events  with  at  least  two  jets  where  the  leading  two 
jets  are  nearly  "back-to-back"  (A(^n  >  150°)  with  ETGet#2)/ET(jet#l)  >  0.8  are 
referred  to  as  "back-to-back"  events  (right) 76 

6-3  Illustration  of  correlations  in  azimuthal  angle  A(^  relative  to  the  direction  of  the 
leading  jet  (highest  Pj  jet)  in  the  event,  jet#l  for  "leading  jet"  events  (left)  and 
"back-to-back"  events  (right)  as  defined  in  Fig.  4.  The  angle  A(|)  =  (j)  -  (fijetttl  is 
the  relative  azimuthal  angle  between  charged  particles  (or  calorimeter  towers) 
and  the  direction  of  jet#l.  On  an  event  by  event  basis,  we  define  "transMAX" 
("transMIN")  to  be  the  maximum  (minimum)  of  the  two  "transverse"  regions, 
60°  <  A(t)<  120°  and  60°  <  -Acj)  <  120°.  "TransMAX"  and  "transMIN"  each  have 
an  area  in  r\-(^  space  of  AT|A(t)  =  47i/6.  The  overall  "transverse"  region  defined  in 
Fig.  3  includes  both  the  "transMAX"  and  the  "transMIN"  region 76 

6-4  Illustration  of  the  topology  of  a  proton-antiproton  collision  in  which  a  "hard" 
parton-parton  collision  has  occurred.  The  "toward"  region  as  defined  in  Fig.  6-1 
contains  the  leading  "jet",  while  the  "away"  region,  on  the  average,  contains  the 


xm 


"away-side"  "jet".  The  "transverse"  region  is  perpendicular  to  tlie  plane  of  the 
hard  2-to-2  scattering  and  is  very  sensitive  to  the  "underlying  event".  For  events 
with  large  initial  or  final-state  radiation  the  "transMAX"  region  defined  in  Fig.6- 
3  would  contain  the  third  jet  while  both  the  "transMAX"  and  "transMIN" 
regions  receive  contributions  from  the  beam-beam  remnants  (see  Fig.  1-1). 
Thus  the  "transMIN"  region  is  very  sensitive  to  the  beam-beam  remnants,  while 
the  "transMAX"  minus  the  "transMIN"  is  very  sensifive  to  initial  and  final-state 
radiation 77 

6-5  Data  at  1 .96  TeV  on  the  density  of  charged  particles,  dNchg/dT|d(|),  with  px  >  0.5 
GeV/c  and  |ri|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  [bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 79 

6-6  Data  at  1 .96  TeV  on  the  density  of  charged  particles,  dNchg/dr|d(t),  with  px  >  0.5 
GeV/c  and  |ti|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Px  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the 
uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the 
particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e. 
generator  level) 80 

6-7  Data  at  1.96  TeV  on  the  PTsum  density  of  charged  particles,  dPTsum/dridc}), 
with  Px  >  0.5  GeV/c  and  |ri|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for 
"leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px 
compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected 
data  (with  statistical  errors  only)  compared  with  the  theory  after  detector 
simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level) 81 

6-8  Data  at  1.96  TeV  on  the  charged  PTsum  density,  dPTsum/dT|d(t),  with  px  >  0.5 
GeV/c  and  |ri|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and 
"transMEN")  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Px  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the 
uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the 
particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e. 
generator  level) 82 
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6-9  Data  at  1 .96  TeV  on  the  average  <pT>  of  charged  particles  with  px  >  0.5  GeV/c 
and  [tiI  <  1  in  the  "transverse"  region  for  "leading  jet"  events  defined  in  Fig.  6-3 
as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA  Tune  A  and 
HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows 
the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the 
statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the 
particle  level  (i.e.  generator  level) 83 

6-10  Data  at  1.96  TeV  on  the  average  maximum  pj,  PTmax,  for  charged  particles 
with  pT  >  0.5  GeV/c  and  JTil  <  1  in  the  "transverse"  region  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 85 

6-11  Data  at  1.96  TeV  on  the  ETsum  density,  dEj/dridcl),  for  particles  with  |ri|  <  1  in 
the  "transMAX"  and  "transMIN"  regions  for  "leading  jet"  events  defined  in  Fig. 
6-3  as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA  Tune  A  and 
HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows 
the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the 
statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the 
particle  level  (i.e.  generator  level) 86 

6-12  Data  at  1.96  TeV  on  the  ETsum  density,  dEx/dridc]),  for  particles  with  |ri|  <  1  in 
the  "transverse"  region  (average  of  "transMAX"  and  "transMIN")  for  "leading 
jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 87 

6-13  Data  at  1.96  TeV  on  the  charged  fraction,  PTsum/ETsum,  in  the  "transverse" 
region  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet 
Pt,  where  PTsum  includes  charged  particles  with  pT  >  0.5  GeV/c  and  hi  <  1  and 
the  ETsum  includes  all  particles  with  |ti|  <  1.  The  data  are  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 88 
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6-14  Data  at  1.96  TeV  on  the  density  of  charged  particles,  dNchg/dr|d(l),  with  pr  >  0.5 
GeV/c  and  |r||  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "back-to- 
back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared 
with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 90 

6-15  Data  at  1.96  TeV  on  the  density  of  charged  particles,  dNchg/dr|d(t),  with  pr  >  0.5 
GeV/c  and  |ri|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Pt  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the 
uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the 
particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e. 
generator  level) 91 

6-16  Data  at  1.96  TeV  on  the  PTsum  density  of  charged  particles,  dPTsum/drid(l), 
with  Pt  >  0.5  GeV/c  and  |ti|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for 
"back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt 
compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected 
data  (with  statistical  errors  only)  compared  with  the  theory  after  detector 
simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level) 92 

6-17  Data  at  1 .96  TeV  on  the  PTsum  density,  dPTsum/dridct),  with  pt  >  0.5  GeV/c  and 
|ri|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and  "transMIN")  for 
"back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt 
compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected 
data  (with  statistical  errors  only)  compared  with  the  theory  after  detector 
simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level) 93 

6-18  Data  at  1.96  TeV  on  the  average  <Pt>  of  charged  particles  with  pt  >  0.5  GeV/c 
and  |ti|  <  1  in  the  "transverse"  region  for  "back-to-back"  events  defined  in  Fig. 
6-3  as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA  Tune  A  and 
HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows 
the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the 
statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the 
particle  level  (i.e.  generator  level) 94 
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6-19  Data  at  1.96  TeV  on  the  average  maximum  pi,  PTmax,  for  charged  particles 
with  pt  >  0.5  GeV/c  and  |t||  <  1  in  the  "transverse"  region  for  "back-to-back" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 95 

6-20  Data  at  1.96  TeV  on  the  ETsum  density,  dEr/dTidct),  for  particles  with  |r||  <  1  in 
the  "transMAX"  and  "transMIN"  regions  for  "back-to-back"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA  Tune  A  and 
HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows 
the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the 
statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the 
particle  level  (i.e.  generator  level) 97 

6-21  Data  at  1.96  TeV  on  the  ETsum  density,  dEx/dTidc]),  for  particles  with  |ri|  <  1  in 
the  "transverse"  region  (average  of  "transMAX"  and  "transMIN")  for  "back-to- 
back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt  compared 
with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 98 

6-22  Data  at  1.96  TeV  on  the  charged  fraction,  PTsum/ETsum,  in  the  "transverse" 
region  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading 
jet  Pt,  where  PTsum  includes  charged  particles  with  pT  >  0.5  GeV/c  and  hi  <  1 
and  the  ETsum  includes  all  particles  with  |ti|  <  1.  The  data  are  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors 
that  include  both  the  statistical  error  and  the  systematic  uncertainty)  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level) 99 

6-23  Data  at  1 .96  TeV  on  the  density  of  charged  particles,  dNchg/dridc]),  with  px  >  0.5 
GeV/c  and  |r||  <  1  in  the  "transMAX"  region  (top),  "transMIN"  region  (middle), 
and  "transverse"  region  (average  of  "transMAX"  and  "transMIN")  (bottom)  for 
"leading  jet"  and  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Pt  compared  with  PYTHIA  Tune  A  and  HERWIG.  The  data  are 
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6-24  Data  at  1.96  TeV  on  charged  PTsum  density  of  charged  particles,  dPTsum/dT|d(t), 
with  Pt  >  0.5  GeV/c  and  |ri|  <  1  in  the  "transMAX"  region  (top),  "transMIN" 
region  (middle),  and  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  (bottom)  for  "leading  jet"  and  "back-to-back"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with  PYTHIA  Tune  A  and 
HER  WIG.  The  data  are  corrected  to  the  particle  level  (with  errors  that  include 
both  the  statistical  error  and  the  systematic  uncertainty)  and  compared  with  the 
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6-25  Figure  6-25:  Data  at  1.96  TeV  on  the  average  <Pt>  of  charged  particles  (top) 
and  the  average  maximum  px,  PTmax,  for  charged  particles  (bottom)  with  pr  > 
0.5  GeV/c  and  |r||  <  1  in  the  "transverse"  region  for  "leading  jet"  and  "back-to- 
back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px  compared 
with  PYTHIA  Tune  A  and  HERWIG.  The  data  are  corrected  to  the  particle 
level  (with  errors  that  include  both  the  statistical  error  and  the  systematic 
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6-26  Data  at  1.96  TeV  on  the  ETsum  density,  dEx/dridcj),  for  particles  with  hi  <  1  in 
the  "transMAX"  region  (top),  "transMIN"  region  (middle),  and  "transverse" 
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6-27  Data  at  1.96  TeV  on  the  charged  fraction,  PTsum/ETsum,  in  the  "transverse" 
region  defined  in  Fig.  6-3  for  "leading  jet"  events  (top)  and  "back-to-back" 
events  (bottom)  as  a  function  of  the  leading  jet  Px,  where  PTsum  includes 
charged  particles  with  px  >  0.5  GeV/c  and  hi  <  1  and  the  ETsum  includes  all 
particles  with  |r||  <  1,  compared  with  PYTHIA  Tune  A  and  HERWIG.  The  data 
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6-28  Data  at  1 .96  TeV  on  the  density  of  charged  particles,  dN/dT|d(t)  and  the  Pxsum 
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7-1  A  comparison  between  the  Run  I  and  Run  II  analysis  on  the  PTsum  density  in 
the  "transverse"  region,  (top)  Data  uncorrected  (with  errors  that  include  both 
statistical  and  the  systematic  uncertainty)  at  1.8  TeV  on  the  charged  PTsum 
density,  dPTsum/dridcl),  with  FT  >  0.5  GeV/c  and  |r||  <  1.0  for  leading  jet  events 
defined  in  Fig.  6-3  as  a  function  of  leading  jet  FT  compared  PYTHIA  set  A, 
PYTHIA  set  B,  and  HER  WIG.  (bottom)  Data  corrected  (with  errors  that  include 
both  statistical  and  the  systematic  uncertainty)  at  1.96  TeV  on  the  charged 
PTsum  density,  dPTsum/dr|d(l),  with  PT  >  0.5  GeV/c  and  |ri|  <  1.0  for  "leading 
jet  events"  defined  in  Fig.  6-3  as  a  function  of  leading  jet  PT  compared  PYTHIA 
Tune  A  and  HERWIG 11 
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CHAPTER  1 
INTRODUCTION 

Elementary  particle  physics  is  concerned  with  the  understanding  of  the 
fundamental  constituents  of  matter.  Our  current  understanding  of  nature  reveals  only  a 
small  number  of  fundamental  particles.  Furthermore,  the  taxonomy  is  greatly  simplified 
because  these  particles  are  perfectly  replicated  in  indistinguishable  copies.  It  should  then 
be  sufficient  to  simply  list  these  particles  and  describe  their  interactions.  Unfortunately, 
due  to  their  size,  we  are  left  to  discover  how  they  interact  via  indirect  measurements. 
Classically  this  requires  the  study  of  decay  rates,  bound  states,  and  scattering 
experiments.  To  fix  the  laws  of  physics  in  a  formulated  phrase,  the  traditional  approach 
is  to  guess  a  form  of  the  interaction  and  compare  the  resulting  theoretical  calculations  to 
experimental  measurements. 

One  important  quantity  measured  by  the  experimentalist  and  calculated  by  the 
theorist  is  the  differential  scattering  cross  section.  Suppose  particles  1  and  2  collide, 
producing  particles  3,4,..., n,  the  cross  section  is  given  by  formula  1.1, 

Y  d'p,  Y  d'p,  ^  (  d'p„ 


da  =  \m\ 
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where  M  is  the  invariant  amplitude  (matrix  element)  for  the  process,  p,  =  [E.  ,  p. )  is  the 

four-momentum  of  particle  /  (mass  mi),  E.  -  yjmf  +  pf  ,  and  S  is  a  statistical  factor 

The  inclusive  jet  cross  section  at  the  Tevatron,  da/dEi  averaged  over  a  small  range 
of  pseudorapidity,  r\  =  -ln(tan  9/2),  is  an  important  object  for  study  because  it  tests 


perturbative  QCD  at  the  highest  Q^  (Q^  =  -q^  =(pi-pf)^)  scale  currently  possible  [1,2]. 
This  measurement  continues  to  substantiate  QCD  and  contributes  to  its  global  data  fitting 
used  to  measure  parton  distribution  functions  [3].  It  is  imperative  that  we  consider  all 
systematic  effects  that  influence  our  interpretation  of  the  measurement. 

One  of  the  contributing  effects  results  from  the  "underlying  event"  that  is  present  in 
all  proton-antiproton  collisions.  Figure  1-1  illustrates  a  "hard"  2-2  parton  scattering  with 
transverse  momentum,  Pjhard-  The  "hard"  scattering  component  is  comprised  of  the  two 
outgoing  jets  plus  initial  and  final-state  radiation.  Removal  of  the  two  outgoing  jets 
leaves  only  the  "underlying  event." 

"Hard"  Scattering/ 

^X  Outgoing  Parton 
/PT(hard) 

Proton       ~^      ~:^       T"^  AntiProton 


Underlying  Event ^^^^^^T^T^T"'  TTnilrrlrinr  Event 


.  Initial-State 
Radiation 


Outgoing  Parton  /      ^ 

Figure  1-1:  Illustration  of  the  way  QCD  Monte-Carlo  models  simulate  a  proton- 
antiproton  collision  in  which  a  "hard"  2-to-2  parton  scattering  with  transverse 
momentum,  p^""' ,  has  occurred.  The  resulting  event  contains  particles  that  originate  from 
the  two  outgoing  partons  (plus  initial  and  final-state  radiation)  and  particles  that  come 
from  the  breakup  of  the  proton  and  antiproton  (i.e.  "beam-beam  remnants").  The 
"underlying  event"  is  everything  except  the  two  outgoing  hard  scattered  "jets"  and 
consists  of  the  "beam-beam  remnants"  plus  initial  and  final-state  radiation.  The  "hard 
scattering"  component  consists  of  the  outgoing  two  jets  plus  initial  and  final-state 
radiation. 

Theoretical  analysis  of  such  a  process  uses  simulations  based  on  QCD  Monte  Carlo 

models.   The  justification  lies  in  the  factorization  theorem  [4-6],  which,  roughly,  states 

that  physical  observables  are  the  product  of  short  distance  functions  and  long  distance 


functions.  The  short  distance  functions  are  calculable  in  perturbation  theory,  where  the 
usual  perturbative  expansion  in  terms  of  Feynman  diagrams  is  used  to  calculate  matrix 
elements.  The  long  distance  functions  are  fit  at  a  scale,  but  their  evolution  to  any  other 
scale  is  also  calculable  in  perturbation  theory.  Of  course,  due  to  the  nature  of  Quantum 
mechanics  it  is  impossible  to  distinguish  initial  state  radiation  from  final  state  radiation. 
However,  in  the  Leading  Log  Approximation  Figure  1-1  can  be  factorized  into  the 
following  subprocesses:  (1)  final  state  emission;  (2)  initial  state  emission;  (3)  the 
elementary  hard  subprocess,  which  can  be  computed  exactly  to  finite  order  in 
perturbation  theory;  and  (4)  The  hadronization  process. 

Refinement  of  each  of  these  subprocesses  extends  the  comparative  scale  of  high 
energy  physics.  The  work  presented  here  extends  to  higher  energies  the  previous 
characterization  of  the  "underlying  event"  [7-13].  Unlike  the  previous  analysis  we 
examine  the  energy  in  the  transverse  regions  and  correct  the  data  back  to  the  particle 
level. 

1.1  Forces  and  Particles 

Four  known  forces  serve  as  the  impetus  of  our  physical  theories:  strong, 
electromagnetic,  weak,  and  gravitational.  Intimately  related  to  these  forces  are  the 
relatively  few  "elementary"  particles  on  which  they  act:  leptons,  quarks,  and  bosons. 

Our  current  understanding  of  gravity  started  with  Sir  Isaac  Newton's  law  of 
universal  gravitation.  This  classical  theory  of  gravity  saw  its  relativistic  generalization  to 
the  general  theory  of  relativity  by  Albert  Einstein.  However,  gravitational  forces  are  the 
weakest,  and  they  are  important  for  massive  bodies  but  negligible  for  nuclear  and 
subnuclear  particles. 


Electrodynamics,  the  theory  describing  electromagnetic  forces,  owes  its  classical 
formulation  to  the  Scottish  physicist  James  Clerk  Maxwell.  It  was  not  until  the  1940's 
that  Tomonaga,  Feynman,  and  Schwinger  would  perfect  the  quantum  theory  of 
electrodynamics.  The  theory  of  weak  interactions  or  flavordynamics  was  originally 
formulated  by  Fermi  in  1933,  but  in  the  1960's  Glashow,  Wienberg,  and  Salam  (GWS) 
put  it  into  its  present  form.  In  the  GWS  model  the  weak  and  electromagnetic  interactions 
are  treated  as  different  manifestations  of  a  single  electroweak  force.  This  electroweak 
force  provides  for  the  attraction  between  charged  particles  and  is  responsible  for  the  [3 
decay  of  nuclei. 

Chromodynamics  emerged  in  the  mid-seventies  to  support  the  pioneering  work  of 
Yukawa  on  the  strong  force.  The  work  of  Gross,  Politzer,  and  Wilczek  on  a  strong  force 
earned  them  recent  recognition  in  the  form  of  the  2004  Nobel  Prize  in  physics.  Strong 
forces  act  only  at  very  small  distances:  they  bind  quarks  into  nucleons  and  nucleons 
together  to  make  nuclei. 

All  of  these  apparently  different  forces  are  each  mediated  by  the  exchange  of  an 
integer-spin  particle,  called  a  boson.  Table  1-1  summarizes  the  properties  of  the  gauge 
bosons  of  the  Standard  Model. 


Table  1-1:  Gauge  bosons  and  forces  of  the  Standard  Model.    There  are  eight  different 
species  of  gluons  each  corresponding  to  a  particular  color  charge. 


Boson 

Force 

Spin 

Charge  [e] 

Mass  [GeV/c^] 

Range  [fm] 

g 

strong 

1 

0 

0 

<~1 

7 

electromagnetic 

1 

0 

0 

00 

W* 

weak 

1 

±1 

80.425±0.038 

=  10-^ 

7^ 

weak 

1 

0 

91.1876±0.0021 

=  10-^ 

The  Standard  Model  [14,  15]  classifies  elementary  particles  as  structureless  at  all 
scales  presently  accessible.  Furthermore,  it  states  that  all  visible  matter  consists  of 
elementary  particles  of  two  kinds:  leptons  and  quarks.  These  particles  are  interpreted  as 
quantum  excitations  of  a  field  and  are  characterized  by  having  spin  Vi  intrinsic  angular 
momentum  in  units  of  h.  They  obey  the  Pauli  Exclusion  Principle  and  are  called 
fermions.  There  are  six  types  of  leptons  and  six  types  of  quarks  which  are  each  grouped 
into  three  generations  according  to  their  mass,  the  properties  of  which  are  summarized  in 
Tables  1-2  and  1-3.  Each  has  an  associated  antiparticle  with  the  same  mass  and  spin  but 
opposite  charge. 


Table  1-2:  Properties  of  leptons. 


Lepton 

Spin 

Charge  [e] 

Mass  [MeV/c^] 

V^  generation 

e' 

Ve 

1/2 
1/2 

-1 
0 

0.51099892±0.00000004 
<3*10-^ 

2"^  generation 

1^" 

1/2 
1/2 

-1 
0 

105.658369±0.000009 
<0.19 

3'^''  generation 

T" 

Vt 

1/2 
1/2 

-1 
0 

1776.99!°o:f, 
<18.2 

Table 

1-3:  Prope 

rties  of  quarks 

Quark 

Spin 

Charge  [e] 

Mass 

f  generation 

u 
d 

1/2 
1/2 

+2/3 
-1/3 

1.5-4  MeV/c' 
4-8  MeV/c^ 

2"''  generation 

c 
s 

1/2 
1/2 

+2/3 
-1/3 

1.15-1.35  GeV/c' 
80-130  MeV/c^ 

3'^''  generation 

t 
b 

1/2 
1/2 

+2/3 
-1/3 

178.0±4.3  GeV/c^ 
4.1-4.4  GeV/c^ 

1.2  The  Standard  Model 

The  Standard  Model  of  strong,  weak,  and  electromagnetic  interactions  is  based  on 
the  symmetry  group  SU(3)cxSU(2)lXU(1)y  [16]  where  the  subscripts  denote  special 
features  of  a  given  symmetry  which  act  on  the  quark  and  lepton  fields.  The  C  in  SU(3)c 
stands  for  color.  Each  quark  has  three  color  components  and  SU(3)c  transforms  them 
into  one  another.  SU(3)c,  the  basis  of  quantum  chromodynamics  (QCD),  is  an  exact 
symmetry  of  nature.  There  are  eight  massless  gluons  which  correspond  to  the  eight 
gauge  fields  of  SU(3)c.  This  description  requires  a  new  quantum  number:  color  charge. 
By  convention,  the  three  colors  are  red,  green,  and  blue  and  each  quark  is  supposed  to 
carry  one  of  these  colors.  The  gluons,  the  quanta  of  color  fields,  also  carry  color.  Quarks 
are  bound  together  in  hadrons  by  the  strong  color  force  via  the  exchange  of  colored 
gluons.  All  observed  hadrons  are  described  in  the  parton  model  as  color  singlet  states 
(referred  to  as  "colorless")  composed  of  three  quarks  (baryons:  qqq)  or  of  a  quark- 
antiquark  pair  (mesons:  qq).  The  quarks  of  these  configurations  are  called  valence 
quarks  because  they  are  responsible  for  the  charge  and  other  quantum  numbers  of 
hadrons. 

The  L  on  SU(2)l,  the  weak  isospin  group,  denotes  the  fact  that  only  left-handed 


components,  ^u  of  spinor  fields,  ^^^ 


V    2    , 


y/ ,  transform  as  doublets  under  that 


group.     The  right-handed  spinor  components   y/,^  = 


V     2     , 


y/   are  isosinglets  under 


SU(2)l;  i.e.,  they  are  unchanged  under  SU(2)l  transformations  and  therefore  do  not 

couple  to  its  three  gauge  fields  which  we  denote  by  W^  =    ''  f-  ^    and  Wj* .   The  fact 

V2 


that  only  left-handed  quarks  and  leptons  couple  to  those  gauge  fields  makes  their  (weak) 
interactions  maximally  parity  violating  [17]. 

The  Y  on  U(1)y  stands  for  weak  hypercharge,  the  charge  associated  with  that 
Abelian  group.    That  gauge  group  has  one  gauge  field  B^  that  couples  to  quarks  and 

leptons  via  their  hypercharge  Y.  It  couples  to  left  and  right  handed  components  of  these 
particles  differently  and  therefore  also  violates  parity,  but  not  maximally.  The 
SU(2)lxU(1)y  part  of  the  standard  model  is  not  an  exact  symmetry.  If  it  were,  the  W^, 

W^ ,  and  B^  would  all  be  massless  gauge  bosons.  That  is  not  the  case.  To  accommodate 

electroweak  phenomenology,  a  scalar  (spin  0)  field  is  introduced  which  breaks  the 
symmetry  SU(2)lxU(1)y  down  to  the  U(l)  symmetry  of  QED.  That  breaking  gives  mass 
to  the  W,f ,  and  the  combination  of  fields  Z,,  =  W^.  cosO^.  -  B,.  sin  0^  which  is  called  the 

Z  boson.  The  orthogonal  combination  A^  =  B^  cos  0^  +  W^  sin  ^„.  remains  massless  and 

is   identified   as   the   photon.      The   angle    0^,   called   the   weak   mixing   angle,   is 

experimentally  found  to  be  sin^^^  =0.23.    That  leads  to  standard  model  predictions 

Massw  =  80GeV/c^  and  Massz  -  91GeV/cl  The  discoveries  of  the  W  [18,  19]  and  Z  [20, 
21]  bosons  by  the  UAl  and  UA2  collaborations  at  Cem  substantiate  the  vahdity  of  the 
Standard  Model.  The  search  for  the  spin  0  Higgs  Boson  continues  to  be  one  of  the  most 
important  problems  to  be  addressed  during  Run  2  at  the  Tevatron  and  future  LHC 
experiments. 

1.3  QCD  and  the  Structure  of  Hadrons 
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QCD  is  a  non-Abelian'  gauge  theory  that  is  based  on  the  SU(3)c  group  of 
transformations  which  relate  quarks  of  different  colors.  The  gauge  bosons  associated 
with  the  eight  group  generators,  known  as  gluons,  can  be  emitted  or  absorbed  by  quarks 
in  transition  in  which  the  color  (but  not  flavor)  can  change.  Since  the  gluons  themselves 
carry  color  they  can  interact  with  each  other  as  well. 

In  quantum  field  theories  like  QED  and  QCD  any  charge  (color  or  electromagnetic) 
is  shielded  by  a  cloud  of  polarized  charges:  a  quark  can  emit  a  gluon  which  can  convert 
into  ^^or  gg  pairs  which  in  turn  can  radiate  gluons  and  we  have  a  branching  tree  of 
quarks  and  gluons  (this  effect  is  called  vacuum  polarization).  Because  of  the  effect  of 
charge  screening  the  charge  one  measures  depends  on  the  distance  (or  wavelength,  or 
transferred  momentum  Q  )  with  which  one  is  probing  the  charge  itself.  We  thus  have  a 
"running"  coupling  constant  which  changes  with  the  transferred  momentum: 

a  =  a(Q')  with    Q"  =-q'- >Q  (1-2) 

where  q  is  the  four-momentum  of  the  virtual  boson  exchanged  between  charges. 

For  both  QED  and  QCD  the  effective  coupling  constant  a  depends  on  the 
momentum  (or  distance)  scale  at  which  it  is  evaluated,  and  takes  the  general  form  of 
equation  1-3. 
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Where  a(0)  for  QED  is  the  fine  structure  constant  and  is  approximately  equal  to  1/137. 
For  QED  it  can  be  shown  that  X(Q^)  takes  the  form  of  equation  1-4. 


'  A  group  or  other  algebraic  object  is  called  non-Abelian  if  the  law  of  commutativity  does  not  always  hold. 
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Here  Nf  is  the  number  of  fundamental  fermions  with  masses  below  '/2|Q|  and  |i  is  the 
mass  of  the  heaviest  fermion  in  the  energy  region  being  considered.  Clearly,  X(Q^)  for 
QED  is  >  0,  and  the  coupling  constant  grows  with  energy.  At  some  energy  scale  the 
coupling  of  QED  becomes  strong  and  perturbation  theory  no  long  applies.  Due  to  this 
behavior,  the  bare  charge  in  QED  is  said  to  be  "ultraviolet  divergent." 

The  coupling  constant  in  QCD  exhibits  the  opposite  behavior.  It  can  be  shown  that 


X(Q  )  for  QCD  takes  the  form  of  equation  1-5. 
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Here  Nf  is  the  number  of  quark  flavors  with  masses  below  V2|Q|  and  |i  is  the  mass  of  the 
heaviest  quark  in  the  energy  region  being  considered,  and  Nc  is  the  number  of  colors.  In 
contrast  to  the  form  of  the  QED  term,  for  6  flavors  and  3  colors  2NrllNc  <  0  and 
therefore  a(Q  )  decreases  with  increasing  momentum  (or  shorter  distances).  Only  in  a 
world  with  more  than  16  quark  flavors  (we  are  safely  below  this  number  at  present 
energies)  is  the  sign  of  X(Q  )  the  same  as  in  QED.  This  results  in  an  antiscreening  of  the 
color  charge:  by  moving  closer  to  the  original  quark  the  amount  of  the  measured  color 
charge  decreases. 

The  QCD  "running  coupling  constant"  (Eq.  1-6)  is  usually  expressed  in  terms  of  a 
parameter,  Aqcd,  that  indicates  the  magnitude  of  the  scale  at  which  as(Q^)  becomes 
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strong;  it  is  determined  experimentally  to  be  about  0.2  GeV. 

Figure  1  -2  compares  the  relative  strength  of  the  strong  and  electromagnetic  forces 
at  different  energy  scales.  With  three  colors  and  six  flavors  we  can  see  that  as(Q  )  in  eq. 
1-6  goes  to  zero  as  Q  goes  to  infinity.  This  results  in  quarks  and  gluons  appearing 
almost  like  free  particles  when  probed  at  very  high  energies  or  short  distances.  This 
behavior  is  called  asymptotic  freedom  and  allows  perturbation  theory  to  be  applied  to 
theoretical  QCD  calculations  to  produce  experimentally  verifiable  predictions  for  hard 
scattering  processes. 

Asymptotic  freedom  is  quickly  overcome  by  the  strong  force  as  color  charges 
separate.  In  contrast,  as  two  electrically  charged  particles  separate,  the  electric  fields 
between  them  quickly  diminish,  allowing  electrons  to  become  unbound  from  nuclei. 
However,  the  gluon  field  lines  associated  with  color  charges  do  not  radially  fan  out  but 
remain  confined  to  a  narrow  cylindrical  region.  This  leads  to  an  interaction  energy  that  is 
proportional  to  the  separation  distance  of  the  sources  of  the  field  lines.  When  quarks 
become  separated,  as  happens  in  high  energy  hadron  collisions,  at  some  point  it  is 
energetically  favorable  for  a  new  quark/antiquark  pair  to  "pop"  out  of  the  vacuum  than  to 
allow  the  quarks  to  separate  further.  As  a  result  of  this,  when  quarks  are  produced  in 
collisions,  instead  of  seeing  the  individual  quarks,  we  see  an  avalanche  of  colorless 
particles  clustered  together  moving  in  roughly  the  same  direction  known  as  jets.  This 
process  is  called  hadronization,  fragmentation  or  string-breaking. 
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Figure  1-2;  Relative  strength  of  the  strong  and  electromagnetic  forces 


1.4  Hadron  Hadron  Interactions 
From  a  phenomenological  point  of  view,  we  can  consider  two  hadrons  colliding  at 
high  energies,  such  as  we  have  at  the  Tevatron  Collider,  to  be  colliding  broad-band 
beams  of  quarks,  antiquaries,  and  gluons.  An  average  hard  scattering  event  consists  of  a 
collection  (or  burst)  of  hadrons  traveling  roughly  in  the  direction  of  each  of  the  initial 
beam  particles  and  two  collections  of  hadrons  with  large  transverse  momentum.  The  two 
large  transverse  momentum  jets  are  roughly  back  to  back  in  azimuthal  angle.  One  can 
use  the  topological  structure  of  hadron-hadron  collisions  to  smdy  the  "underlying  event" 
[13,  22-24].  Of  course,  from  a  certain  point  of  view  there  is  no  such  thing  as  an 
"underlying  event"  in  a  proton-antiproton  collision.  There  is  only  an  "event"  and  one 
cannot  say  where  a  given  particle  in  the  event  originated.  On  the  other  hand,  hard 
scattering  collider  "jet"  events  have  a  distinct  topology.  On  the  average,  the  outgoing 
hadrons  "remember"  the  underlying  2-to-2  hard  scattering  subprocess. 
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1 .4. 1  Parton  Model  and  Large  Px  Processes 
Hadronic  collisions  which  involve  a  short  distance  scattering  can  be  described,  in 
first  approximation,  by  the  Parton  Model.  Since  the  transferred  momentum  is  the 
conjugate  variable  of  the  (quark-quark)  impact  parameter,  larger  Q  implies  that  partons 
have  scattered  at  small  distances  where  tts  is  small.  In  order  to  be  able  to  apply 
perturbation  theory  one  needs  a  momentum  transfer  of  about  10  GeV:  from  the 
uncertainty  principle  we  calculate  that  the  associated  distance  is  in  the  order  of  lO'^m. 
From  the  experimental  point  of  view  there  is  no  certain  way  to  separate  such  rare  events 
but  intuition  suggests  that  large  pt  final  state  particles  should  be  a  good  indication  and 
experiments  fully  support  this. 

We  can  picture  this  scattering  process  as  a  sequence  of  three  different  phases 
occurring  at  different  time  scales.  Partons  approach  each  other  carrying  a  fraction,  x,  of 
the  momentum  of  their  parent  hadrons.  The  transverse  momenta  of  the  partons  are 
neglected.  "Parton  Distribution  Functions"  (PDF)  f\x,/i^)  are  so  introduced  giving  the 

probability  for  parton  /  to  have  fractional  momentum  between  x  and  x  +  etc  (|j.  is  a 
factorization  scale).  Figure  1-3  is  an  example  of  such  parton  distribution  functions. 
These  functions  are  extracted  from  experimental  data  in  deep  inelastic  scattering  (DIS)  of 
leptons  on  nucleons  [25].  These  parton  distribution  functions  of  the  proton  are  a  result  of 
the  work  done  by  the  HI  [26-28]  and  ZEUS  [29]  collaborations  at  HERA,  and  of  the 
inclusive  jet  distribution  by  DO  [30]  and  CDF  [31]  collaborations  at  the  Tevatron.  As  the 
data  only  cover  a  finite  range  of  Q^,  the  evolution  of  such  functions  with  Q^  is  computed, 
using  perturbation  theory,  with  the  Altarelli-Parisi  equation  [32]. 
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Figure  1-3:  The  parton  structure  functions  extracted  from  an  analysis  of  deep  inelastic 
scattering  data  at  Q^=10GeVl 


A  hard  collision  then  takes  place  between  a  pair  of  partons  regarded  as  free 
particles.  Predictions  for  jet  production  are  given  by  folding  the  parton  distribution 
functions  with  perturbatively  calculated  "two-body"  scattering  cross  sections  &.j .   Any 


cross  section  of  interest  is  calculated  using  equation  (1-7). 
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Figure  1-4  shows  the  representation  of  an  elementary  "two-body"  interaction 
between  two  partons  in  a  pp  collision  producing  a  di-jet  event. 
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Figure  1-4:    Hard  "two-body"  parton  interaction  producing  a  di-jet  event  in  a  proton- 
antiproton  collision. 

The  factorization  scale  |i  discriminates  whether  a  parton,  inside  the  incoming 
hadron,  takes  part  or  not  in  the  hard  scattering:  if  the  momentum  of  a  parton  is  greater 
than  the  scale  \i,  it  contributes  to  the  short-distance  cross  section  (as  the  partons  /  and  7  in 
fig.  1-4);  if  its  momentum  is  less  than  the  scale  |i,  it  is  considered  part  of  the  hadron 
structure  not  involved  in  the  hard  interaction  (spectator  parton). 

1.4.2  The  "Underlying  Event"  in  Proton-Antiproton  Collisions:  Pythia  and  Herwig 

Fig.  1-1  illustrates  the  way  QCD  Monte-Carlo  models  simulate  a  proton-antiproton 
collision  in  which  a  "hard"  2-to-2  parton  scattering  with  transverse  momentum,  PT(hard), 
has  occurred.  The  resulting  event  contains  particles  that  originate  from  the  two  outgoing 
partons  (plus  initial  and  final-state  radiation)  and  particles  that  come  from  the  breakup  of 
the  proton  and  antiproton  (i.e.,  "beam-beam  remnants").  The  "underlying  event"  is 
everything  except  the  two  outgoing  hard  scattered  "jets"  and  receives  contributions  from 
the  "beam-beam  remnants"  plus  initial  and  final-state  radiation.  The  "hard  scattering" 
component  consists  of  the  outgoing  two  jets  plus  initial  and  final-state  radiation.  Any 
measured  observable  of  the  underlying  event  necessarily  receives  contributions  from 
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initial  and  final  state  radiation.    It  is  possible  to  reduce  these  contributions  by  placing 
constraints  on  event  topology  and  jet  energy. 

The  "beam-beam  remnants"  are  what  is  left  over  after  a  parton  is  knocked  out  of 
each  of  the  initial  two  beam  hadrons.  It  is  the  reason  hadron-hadron  collisions  are  more 
"messy"  than  electron-positron  annihilations  and  no  one  really  knows  how  it  should  be 
modeled.  For  the  QCD  Monte-Carlo  models  the  "beam-beam  remnants"  are  an  important 
component  of  the  "underlying  event."  Also,  it  is  possible  that  multiple  parton  scattering 
contributes  to  the  "underiying  event."  Figure  1-5  shows  the  way  PYTHIA  [33]  models 
the  "underlying  event"  in  proton-antiproton  collision  by  including  multiple  parton 
interactions.  In  addition  to  the  hard  2-to-2  parton-parton  scattering  and  the  "beam-beam 
remnants,"  sometimes  there  is  a  second  "semi-hard"  2-to-2  parton-parton  scattering  that 
contributes  particles  to  the  "underlying  event." 

For  the  hadronization  process  we  have  three  types  of  non-perturbative  contributions 
to  consider:  (1)  representation  of  the  incoming  partons  as  constituents  of  the  incident 
hadrons  via  parton  distribution  functions;  (2)  the  conversion  of  the  emitted  partons  into 
outgoing  hadrons  using  quark  and  gluon  fragmentation  functions;  and  (3)  the  "soft" 
component  to  the  "underlying  event"  generated  by  spectator  partons.  The  "underlying 
event"  receives  contributions  from  the  original   pp   system  ("beam-beam  remnants"), 

initial  and  final-state  radiation,  and  possibly  hadrons  resulting  from  multiple  parton 
interactions,  as  in  Figure  1-5. 
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Fig.  1-5.  Illustration  of  the  way  PYTHIA  models  the  "underlying  event"  in  proton- 
antiproton  collision  by  including  multiple  parton  interactions.  In  addition  to  the  hard  2-to- 
2  parton-parton  scattering  with  transverse  momentum,  PxChard),  there  is  a  second  "semi- 
hard" 2-to-2  parton-parton  scattering  that  contributes  particles  to  the  "underlying  event". 

Both  of  the  QCD  Monte  Carlo  models,  HERWIG  [34]  and  PYTHIA,  include  a 

"soft"  underlying  event  structure  that  is  modeled  by  a  parameterization  of  minimum  bias 

data,  and  creates  jets  from  the  beam  remnants.     This  is  an  imperfect  model  of  the 

"underlying  event,"  since  it  always  contains  particles,  yet  the  soft  underlying  event  can 

be  absent,  as  shown  by  the  finite  survival  probabilities  for  inelastic  events  with  large 

rapidity  gaps  [35-38].     PYHTIA  has  added  multiple  parton  interactions  to  enhance  the 

activity  of  the  "underiying  event,"  and  Tune  A  [22-23]  was  specifically  tuned  to  fit  the 

"underlying  event"  in  the  Run  I  data. 

The  QCD  perturbative  2-to-2  parton-parton  differential  cross  section  diverges  as  the 
transverse  momentum  of  the  scattering,  Pj,  goes  to  zero  (see  figure  1-1).  PYTHIA  uses  a 
tunable  parameter  to  prevent  divergences  at  low  Pj.  Tune  A  was  tuned  to  fit  the  Run  I 
data,  and  the  low  Pt  region  was  an  area  of  focus.  HERWIG  does  not  allow  for  this,  and  a 
suitable  Pt  cutoff  most  be  chosen.  We  use  Pt  >  5  GeV/c  for  all  2-to2  hard  scattering 
events  in  HERWIG. 


CHAPTER  2 
ACCELERATOR  AND  DETECTOR 

The  Tevatron  [39]  pp  Collider  is  currently  the  world's  highest  energy  particle 
accelerator  in  operation.  It  is  the  largest  in  a  chain  of  five  accelerators  at  the  Fermi 
National  Accelerator  Laboratory  (FNAL,  Fermilab)  and  is  capable  of  producing  proton- 
antiproton  collisions  at  a  center  of  mass  energy  V^  =1.96  TeV.  The  Collider  Detector  at 
Fermilab  (CDF)  [40]  and  DO  [41]  are  the  multipurpose  detectors  built  at  collision  points 
to  exploit  physics  at  the  Tevatron.  The  analysis  presented  in  this  dissertation  is  based  on 
the  data  sample  collected  by  CDF  during  the  2001-2004  (Run  II)  running  period  of  the 
Tevatron. 

2. 1  The  Accelerator  Complex 

The  Accelerator  Complex  at  FNAL  (see  Figure  2-1)  uses  multiple  stages  of 
acceleration  to  achieve  proton-antiproton  collisions  at  a  center  of  mass  energy  Vi^= 1.96 
TeV.  The  protons  used  in  the  collisions  originate  from  ionized  Hydrogen  gas  molecules. 
These  H"  ions  are  first  accelerated  to  750  KeV  in  the  Cockcroft-Walton  accelerator.  They 
are  then  fed  into  the  Linac,  a  150m  long  series  of  nine  radio-frequency  (RF)  cavities 
which  produce  an  electric  field  that  rapidly  changes  direction.  In  this  linear  accelerator 
the  H"  ions  are  brought  to  400  MeV.  Subsequent  to  this  stage,  the  beam  is  focused  and 
made  to  collide  with  a  thin  fixed  carbon  target  which  affects  the  loss  of  two  electrons  per 
ion.  The  denuded  H"  ions  are  now  the  protons  that  will  ultimately  be  collided  or  used  to 
make  the  anti-protons  at  the  Target  Station. 
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The  75.5  m  radius  Booster  is  a  fast  cycling  proton  synchrotron  of  conventional 
magnets  (used  to  steer  and  focus  the  beam)  and  an  RF  cavity  (used  to  accelerate  the 
beam).  The  accelerated  protons  leave  with  kinetic  energy  of  8  GeV  and  are  then  injected 
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TEVATRON 
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Figure  2-1 :  Overview  of  the  accelerator  complex  at  Fermilab.  H-  ions  are  injected  into 
the  linac  from  the  Cockcroft-Walton,  to  travel  to  the  Booster,  then  to  the  Main  Ring,  and 
finally  to  the  Tevatron.  Some  protons  are  extracted  from  the  Main  Ring  and  are  used  to 
make  anti-protons.  The  anti-protons  are  re-injected  into  the  Main  Ring  and  then  into  the 
Tevatron.  The  final  pp  center  of  mass  energy  is  V^=1.96  TeV. 

into  the  Main  Ring.  These  protons  are  then  further  accelerated  to  150  GeV  in  the  Main 
Injector,  and  finally  brought  to  980  GeV  by  the  Tevatron.  The  Tevatron  was  the  world's 
first  superconducting  synchrotron.  The  beam  is  guided  around  the  closed  path  by  dipole 
magnets.  As  the  beam  energy  is  ramped  up  by  RF  cavities  from  150  GeV  to  980  GeV, 
the  bending  magnetic  fields  and  the  RF  frequency  must  be  carefully  synchronized  to 
ensure  beam  stability.  The  transverse  motion  of  the  beam  is  stabilized  by  quadrupole 
magnets  that  take  advantage  of  magnetic  field  gradient  technology. 
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Anti-protons  are  produced  by  sending  120  GeV  protons  from  the  Main  Injector  to  a 
nickel  target.  From  the  resulting  shower  of  particles,  antiprotons  of  around  8GeV  are 
selected  and  sent  to  the  Debuncher  and  Accumulator  Rings  where  RF  and  stochastic 
cooling  systems  are  used  in  the  momentum  stacking  process.  Once  a  'stack'  has  been 
collected  the  antiprotons  are  sent  back  to  the  Main  Injector,  accelerated  to  150  GeV  and 
put  into  the  Tevatron,  circulating  counter  to  the  proton  bunches.  For  every  million 
protons  that  hit  the  target,  only  about  twenty  8  GeV  anti-protons  will  be  stacked  into  the 
Accumulator. 

The  Recycler  is  placed  directly  above  the  Main  Injector  beamline.  This  dual 
function  system  serves  as  a  post  Accumulator  storage  ring  and  as  a  recycler  for  the 
antiprotons  left  over  at  the  end  of  a  store'.  These  recycled  anti-protons  can  be  mixed 
with  those  from  the  Accumulator  and  then  accelerated  to  150  GeV  in  the  Main  Injector 
and  then  injected  into  the  Tevatron. 

Once  both  beams  are  at  the  maximum  energy  they  are  focused  and  brought  to 
collision  at  the  two  interaction  points,  one  of  which  is  at  the  center  of  the  CDF  detector. 
The  luminous  region  has  a  Gaussian  dispersion  of  around  30|i,m  transverse  to  the  beam 
direction,  and  a  length  along  the  beam  direction  of  around  30cm. 

The  beams  typically  circulate  for  12-18  hours  during  which  time  the  luminosity 
falls  approximately  an  order  of  magnitude.  During  this  time  antiproton  are  continuously 
stacked.  When  the  stack  is  sufficiently  large  and  the  luminosity  has  significantly 
decayed,  the  beam  in  the  Tevatron  is  dumped  and  new  bunches  are  injected. 


A  store  refers  to  the  period  of  time  when  proton-antiproton  collisions  are  taking  place.  Stores  can  be 
selectively  terminated  when  luminosity  has  reached  a  minimum  or  a  new  stack  of  antiprotons  is  ready. 
Altematively,  a  store  can  end  with  a  magnet  quench  or  other  problem. 
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The  instantaneous  luminosity  L  is  given  by  equation  2-1 

fN„N- 


(2-1) 


where  /  is  the  frequency  of  bunch  crossings,  Np  and  A^- ,  are  the  number  of  protons  and 

antiprotons,  respectively,  per  bunch,  and  A  is  the  effective  area  of  the  crossing  beams. 
The  current  status  of  the  luminosity  is  shown  in  figure  2-2,  and  the  integrated  luminosity 
delivered  to  tape  is  shown  in  Fig.  2-3. 

Parts  of  the  Fermilab  accelerator  complex  are  20  years  old  and  there  have  been 
some  setbacks  since  the  upgrades  to  take  the  Tevatron  from  1.8  TeV  to  1.96  TeV. 
However,  the  initial  instantaneous  luminosity  of  stores  has  been  steadily  increasing.  The 
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Figure  2-2:  Run  11  instantaneous  initial  luminosity 
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Figure  2-3:  Run  II  integrated  luminosity 

projected  goal  for  all  of  Run  II  is  4-8  fb"'  by  the  time  LHC  is  ready  to  begin  taking  data 
(circa  2009). 

The  total  number  of  events  n  in  a  scattering  process  is  proportional  to  the 
luminosity  and  the  cross  section  a  of  the  process. 


n  =  La 


2-2 


We  can  get  a  rough  sense  of  the  reach  for  new  physics  and  the  challenge  of  enhancing 
signal  and  suppressing  background  by  considering  the  following  examples.  At  a  center- 
of-energy  of  1 .96  TeV,  we  have: 

CJypp  -^  anything)-  15mb  2-3 

<j{pp  -^  tt  +  anything)  ~  6pb  2-4 

With  about  Ifb"'  of  delivered  luminosity  we  should  have  seen  6000  top  events.  However, 

not  every  second  of  delivered  luminosity  is  observed.  Moreover,  due  to  finite  capabilities 
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in  data  storage  not  every  observed  event  can  be  recorded.    The  task  of  observing  the 
events  fall  on  the  Detector  and  the  job  of  selecting  events  to  be  recorded  is  given  to  the 


trigger  system. 


2.2  The  Collider  Detector  at  Fermilab 


The  Collider  Detector  at  Fermilab  (CDF)  is  a  large,  multilayered  general  purpose 
detector  designed  to  study  a  wide  range  of  processes  occurring  in  proton-antiproton 
collisions.  Figure  2-4  is  a  schematic  drawing  of  the  approximately  5000  ton,  10  m  high, 
and  27  m  long  detector.  The  CDF  is  cylindrically  symmetric  about  the  beam  axis  and 
has  a  forward-backward  symmetry  in  its  tracking,  calorimetric,  and  muon  systems. 
Figure  2-5  shows  an  elevated  view  in  which  the  tracking  system  is  seen  to  be  contained  in 
a  solenoid  coil.  The  calorimetry  and  muon  systems  are  outside  the  solenoid.  These  sub- 
systems are  described  in  more  detail  below. 


Figure  2-4:  Solid  cutaway  view  of  the  CDF  H  detector. 
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Figure  2-5:  Elevation  view  of  the  CDF  II  detector. 
2.2.1  The  CDF  Coordinate  System 
The  geometric  center  of  the  detector  serves  as  the  nominal  interaction  point  (0,0,0). 
Figure  2-6  shows  the  overall  CDF  coordinate  frame,  in  both  cylindrical  and  Rectangular 
systems.  The  positive  z-axis  corresponds  to  the  proton  beam  direction,  the  positive  y- 
axis  points  vertically  upward,  and  the  positive  x-axis  points  radially  outward  from  the 
center  of  the  Tevatron  ring.  In  the  cylindrical  system  the  azimuthal  angle  ^  is  measured 
about  the  beam  axis  from  the  positive  x-axis.  The  polar  angle  9  is  defined  as  the  angle 
measured  from  the  positive  z-axis. 
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Figure  2-6:  The  CDF  coordinate  system. 

The  incident  proton  and  anti-proton  have  no  transverse  momentum  and  they  have 
equal  and  opposite  longitudinal  momentum,  therefore  the  total  momentum  of  the 
products  of  the  collision  would  sum  to  zero  in  a  full  4n  solid  angle  experiment. 
However,  some  space  must  be  left  for  the  beam  pipe.  Those  particles  from  the  collision 
that  travel  at  very  small  angles,  two  degrees  or  less  (such  as  hadrons  from  spectator  quark 
hadronization),  will  fly  down  the  beam  pipe  completely  missing  the  detectors.  Such 
unmeasured  particles  will  not  carry  much  transverse  momentum,  but  they  may  carry 
significant  amounts  of  longitudinal  momentum.  For  this  reason,  longitudinal  momentum 
will  not  balance  in  the  detector,  but  the  transverse  momentum  will,  to  the  detectors 
accuracy. 

For  these  reasons,  rather  than  using  the  total  energy  E  and  total  momentum  p,  we 
generally     use     the     transverse     energy      £,  =£xsin(^)      and     the     transverse 

momentum p,  =  px sin(^) .  In  the  large  energy  collisions  found  at  CDF,  the  Et  and  pt  of 
particles  in  the  event  are  nearly  equally.  However,  by  convention  we  use  Et  when 
referencing  energy  deposited  in  the  calorimeters,  so  that  it  is  understood  that  the  angle  6 
refers  to  the  geometric  center  of  the  detector.  Alternatively,  pt  usually  refers  to  the 
transverse  momentum  of  a  particle  determined  in  a  tracking  chamber  so  that  the  angle  0 
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reflects  the  true  interaction  vertex.  It  is  common  to  use  the  pseudorapidity  in  place  of  the 
polar  angle  0. 

The  "natural"  kinematic  variables  for  hadron  coUisions  are  pseudorapidity, 
transverse  momentum,  and  azimuthal  angle  since  the  shapes  of  their  distributions  are 
invariant  under  a  Lorentz  boost.  Transverse  momentum  and  azimuthal  angle  are 
invariant  to  Loretaz  transformations  along  the  z  axis  and  the  pseudorapidity  is  simply 
additive. 

The  major  components  of  the  CDF  detector  are  arranged  cylindrically  around  the 
interaction  point.  Closet  to  the  beamline  are  the  layers  of  silicon,  providing  high- 
precision  tracking  and  vertexing  in  the  r-^  plane.  Next  is  a  wire  drift  chamber  that 
provides  measurements  of  momentum  and  spatial  parameters  of  a  particles  trajectory 
(track).  The  tracking  subsystem  is  embedded  inside  a  superconducting  solenoidal  magnet 
that  produces  a  1.41  Tesla  magnetic  field.  Energy  measurements  of  jets,  electrons, 
photons,  and  hadrons  are  made  by  the  combined  calorimetry  systems:  central,  plug  and 
forward  electromagnetic  (EM)  and  hadronic  calorimeters  (HAD).  Muons  are  identified 
by  the  presence  of  a  track  in  the  muon  chambers  matched  to  a  track  in  the  central  tracking 
chamber.  Because  of  the  long  lifetime  and  high  penetration  of  muons,  the  muon 
chambers  are  placed  outside  of  the  hadronic  calorimeters,  after  a  steel  absorber  to 
eliminate  any  electromagnetic  and  hadronic  showers.  In  front  of  the  backward  and 
forward  calorimeters  is  a  plane  of  scintillation  counters  called  "Beam-Beam  Counters" 
(BBC).  They  provide  a  minimum  bias  trigger  for  the  detector  and  are  also  used  as  the 
primary  luminosity  monitor.  The  layout  of  these  detectors  is  shown  in  figure  2-7,  which 
depicts  one  quadrant  of  the  cross-section  through  the  detector. 
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Figure  2-7:  A  quarter  of  the  CDF  detector.  Only  the  central  and  end-plug  subsystems  are 
shown. 

2.2.2  Tracking 

The  CDF  uses  silicon  strip  detectors  (LOO  +  SVXII  -i-  ISL)  and  a  drift  chamber 

(COT)  for  charged-particle  track  reconstruction  and  vertex  finding.  The  tracking  systems 

are  inside  a  superconducting  solenoid  of  radius  1.5m  that  provides  a  1.41  Tesla  magnetic 

field  parallel  to  the  beam  axis.  The  magnetic  flux  is  retumed  through  a  steel  yoke.  The 

tracking  volume  and  the  endplug  calorimeters  are  shown  in  figure  2-8.    The  yoke  also 

functions  as  a  support  to  the  calorimeters  located  radially  outside  the  solenoid.  The 

silicon  system  and  drift  chambers  were  redesigned  and  completely  rebuilt  between  Run  I 

and  Run  EI  of  the  Tevatron. 
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Figure  2-8:  The  CDF  11  tracking  volume. 

In  the  silicon  tracker,  as  charged  particles  move  through  the  'depletion  layer' 
created  in  a  biased  p-n  semiconductor  junction,  they  create  electron-hole  pairs  that  drift  to 
be  collected  at  the  surfaces.  This  induces  a  signal  on  metal  stripes  that  have  been 
deposited  on  the  surface  and  connected  to  readout  amplifiers.  Figure  2-9  shows  a 
schematic  layout  of  the  silicon  tracking  system. 

Layer  00  (LOO)  is  mounted  on  the  beam  pipe,  1.6  cm  from  the  beam  axis,  and 
consists  of  8  single  sided  microstrip  silicon  detectors.  These  detectors  cover  the  beam 
pipe  for  about  40  cm  in  each  direction  about  z  =  0  [42].  Outside  the  LOO,  the  'Silicon 
Vertex  Detector'  (SVXII)  occupies  the  volume  between  2.4  and  10.6  cm  from  the  beam 
axis  and  covers  a  total  length  of  96  cm  along  the  z  coordinate.  The  SVXII  system 
consists  of  five  double  sided  microstrip  silicon  layers  [40].  Three  of  these  layers  provide 
track  position  in  the  r-(\)  plane  from  the  readout  of  one  side  (microstrips  parallel  to  the 
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beam  axis)  while  the  z  coordinate  is  determined  by  the  other  side  (microstrip 
perpendicular  to  the  beam  axis).  The  other  two  other  layers  have  their  microstrips  tilted  in 
such  a  way  as  to  provide  a  3-D  track  reconstruction  with  an  approximately  uniform 
efficiency.  The  Intermediate  Silicon  Layers  (ISL)  Detector  is  placed  in  the  region 
between  SVXII  and  the  central  tracking  system  [40].  Figure  2-9  shows  the  position  of 
the  ISL  detector.  The  layers  are  in  the  radial  range  20  <  r  <  30  cm  and  extend  to  |z|  =  65 
cm  for  the  inner  layer  and  |z|  =  87.5  cm  for  the  outer  layer.  The  ISL  covers  the  range  of 
pseudorapidity  for  It)!  ^  2.  Figure  2-10  shows  an  end  view  of  the  three  components  of  the 
silicon  microstrip  detector  system. 
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Figure  2-9:     Schematic  layout  of  the  silicon  tracking  system.    The  innermost  layer, 
LayerOO  consists  of  6  sensors  in  z. 
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The  Central  Outer  Tracker  (COT)  is  a  large  cylindrical  drift  chamber  used  to 
determine  precise  position  measurements  [40].  The  tracking  at  large  radii  in  the  central 
rapidity  region  (|'n|<l)  is  done  with  a  large  open  cell,  cylindrical  drift  chamber  using  a 
readout  that  can  record  multiple  hits  from  each  sense  wire.  The  active  volume  of  the 
COT  spans  310  cm  in  the  beam  (axial)  direction,  z;  between  43.4  cm  and  132.3  cm  in 
radius,  r;  and  the  entire  azimuth,  ^.  The  COT  provides  96  measurement  layers,  organized 
into  alternating  axial  and  ±2°  stereo  superlayers.  Sense  wires  and  potential  wires  are 
alternated  and  arranged  in  8  'superlayers'  as  shown  in  Figure  2-11,  each  consisting  of  12 
layers  of  sense  wires.  Within  each  superlayer  are  'cells',  bounded  by  field-shaping 
sheets.  The  cells  are  angled  at  35°  to  the  radial  direction  to  compensate  for  the  Lorentz 
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Figure  2-10:  End  view  of  the  three  components  of  the  silicon  microstrip  detector  system. 
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Figure  2-11:  The  COT  sense  wires  and  potential  wires  are  alternated  and  arranged  in  8 
'superlayers'. 

angle  of  the  drifting  charged  particles.  There  is  a  'spacer'  at  z  =  0  that  results  in  a  lower 
tracking  trigger  efficiency  at  T)  =  0.  The  chamber  is  filled  with  a  50:50  mixture  of 
Argon/Ethane  and  a  small  amount  of  alcohol.  The  hit  position  resolution  is 
approximately  140|im  and  the  momentum  resolution cr(/7j^)//?^  =  0.0015(Gey/c)~' .  A 
reconstructed  track  provides  accurate  information  in  the  r-^  view  for  measurement  of 
transverse  momentum,  pj,  and  substantially  less  accurate  information  in  the  r-z  view  for 
the  measurement  of  r|. 

2.2.3  Calorimeters 

As  charged  particles  progress  through  the  calorimeters  they  interact  and  develop 

characteristic  'showers'.    Different  size  and  thickness  plastic  scintillator  and  absorber 

layers  are  alternatively  stacked  forming  the  electromagnetic  calorimeters  (allowing  for 

the  energy  measurement  of  photons  and  electrons)  and  the  hadronic  calorimeters 
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(measuring  hadron  energies).  The  primary  particle  produces  a  shower  of  secondary 
particles  inside  the  absorber.  The  shower  particles  deposit  a  fraction  of  their  energy  in 
the  sampling  material  producing  a  light  signal  read  by  photomultipliers  (PMTs)  through 
wavelength  shifting  (WLS)  light  guides  or  optical  fibers. 

As  can  be  seen  in  figure  2-7,  the  CDF  calorimeters  are  physically  separated  into  two 
sections:  the  central  region,  cylindrical  about  the  beam  line  and  covering  |r||  <  1;  and  the 
forward  or  'end  plug'  regions,  covering  1.1  <  [tiI  <  3.6.  The  principal  components  of  the 
central  calorimeter  are  the  central  electromagnetic  (CEM)  [43]  and  the  central  hadronic 
(CHA)  [44]  compartments.  Both  the  CEM  and  CHA  are  retained  from  Run  I.  They  are 
segmented  in  r]  and  (])  with  a  projective  "tower"  geometry,  shown  in  action  in  Fig.  2-12. 
In  each  tower  the  electromagnetic  compartment  is  backed  by  the  hadronic  one,  both 
readout  by  different  PMTs. 

The  central  calorimeter  is  azimuthally  arranged  in  48  physically  separated  15°  wide 
wedges  each  segmented  in  T]  into  ten  towers,  shown  in  figure  2-13.  Each  tower  subtends 
0.1 1  X  15°  in  (j)  X  T).  The  central  electromagnetic  calorimeter  (CEM)  is  overlapped  by  a 
hadronic  section  split  into  two  parts,  the  central  hadronic  calorimeter  (CHA)  and  the  wall 
hadronic  calorimeter  (WHA).  The  CEM  covers  0  <  |ri|  <  1.1  and  uses  lead  sheets 
interspersed  with  polysterene  scintillator  as  the  active  medium  and  employs  phototube 
readout.  The  CHA  covers  0  <  |ti|  <  0.9  and  the  WHA  covers  0  <  |ti|  <  1.3.  Both  hadronic 
calorimeters  use  steel  absorbers  interspersed  with  acrylic  scintillator  as  the  active 
medium. 

Located  six  radiation  lengths  deep  in  the  CEM  calorimeters,  corresponding  to  the 
depth  at  which  showers  typically  reach  their  maximum  transverse  extent  is  the  Central 
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Electromagnetic  Strip  Detector  (CES).  The  CES  uses  proportional  strip  and  wire 
counters  in  a  fine-grained  array,  as  shown  in  Figure  2-14,  to  provide  precise  position 
(about  2  mm  resolution)  and  shape  information  for  electromagnetic  cascades. 


Figure  2-12:  Calorimeter  tower  segmentation  in  r\-^  space. 

A  further  component  of  the  central  calorimeters  is  the  central  pre-radiator  (CPR),  a 
set  of  proportional  chambers  between  the  CEM  and  the  magnet  designed  to  help  separate 
electrons  and  pions. 

The  plug  calorimeters  consist  of  the  plug  electromagnetic  calorimeter  (PEM)  [45], 
newly  built  for  the  CDF  Run  11,  and  the  plug  hadronic  calorimeter  (PHA).  Like  the 
CEM,  the  PEM  consists  of  a  stack  of  lead  and  scintillator  sheets  read  out  by  phototubes. 
At  6  times  the  radiation  length  in  the  PEM  is  the  plug  shower  maximum  detector 
(PES)[46]. 
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Figure  2-13:  CEM/CES/CHA  wedge. 
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Figure  2-14:  CES  strip  and  wire. 
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Finally,  the  first  layer  of  the  PEM  is  read  out  separately  and  referred  to  as  the  plug 
pre-radiator  (PPR).  The  PPR  can  help  to  distinguish  between  electron/photons  and 
hadrons  by  indicating  the  extent  to  which  the  particle  shower  has  already  developed  at  the 
face  of  the  calorimeter. 

2.3  The  CDF  Trigger  System 

The  inelastic  proton-antiproton  cross  section  at  V^=  1.96  TeV  is  about  60  mb 
(-60x10'  cm  ).  For  a  typical  instantaneous  luminosity  of  about  1.2x10  cm"  /s  we  get 
~7.6  million  inelastic  collisions  per  second  at  CDF.  The  CDF  readout  electronics  and 
event  storage  system  is  not  capable  of  recording  events  at  such  a  high  rate.  Moreover, 
most  of  these  events  do  not  present  a  significant  interest  for  the  CDF  physics  program. 
The  trigger  system  is  used  to  select  an  event  rate  of  75  Hz  from  the  7.6  MHz  (132  ns 
crossing)  beam  crossing  rate  [40].  The  event  rate  is  such  that  it  is  necessary  to  filter 
physically  interesting  events  to  be  written  to  tape,  and  this  achieved  through  a  three-level 
trigger  system,  designed  to  be  'deadtimeless'. 

The  Level- 1  trigger  is  achieved  with  hardware.  Based  on  preliminary  information 
from  tracking,  calorimetry,  and  muon  systems,  the  output  of  the  7.6  MHz  Synchronous 
pipeline  with  a  5544ns  latency  at  the  first  level  of  the  trigger  is  used  to  limit  the  rate  of 
accepted  events  to  <50  kHz.  Each  next  trigger  level  examines  fewer  events  but  in  greater 
detail.  At  the  next  trigger,  with  more  refined  information  and  additional  tracking 
information  from  the  silicon  detector,  the  -20  |is  latency,  asynchronous  2  stage  pipeline, 
reduces  the  acceptance  further  to  -  300  Hz.  The  Level  2  algorithm  uses  the  information 
about  high  momentum  tracks  and  clustered  calorimeter  energy.  If  the  accept  decision  is 
made  by  the  trigger,  then  the  information  form  all  subsystems  is  read  out  and  passed  on 
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to  Level  3.  At  level-2  many  triggers  are  prescaled  to  reduce  the  total  acceptance  rates  to 
a  maximum  of  about  20  Hz  that  is  the  maximum  that  Level-3  can  handle.  This  means 
that  a  predefined  fraction  of  events  that  passed  the  trigger  are  considered  to  fail  it. 
Prescaling  is  sometimes  preferable  over  making  tighter  trigger  cuts.  This  method  allows 
us  to  record  as  many  rare  events  as  possible  while  still  accepting  other  data  at  a 
reasonable  rate.  Some  triggers  can  be  dynamically  prescaled.  A  dynamical  prescale  can 
be  changed  during  the  course  of  the  run  depending  on  the  instantaneous  luminosity:  it 
will  be  large  when  the  luminosity  is  high  and  small  when  the  luminosity  is  low. 

Level-3  consists  of  the  event  builder  (EVB)  and  the  Level-3  farm.  The  EVB 
assembles  event  fragments  from  level- 1  and  level-2  into  complete  events,  and  then  the 
Level-3  farm  runs  a  version  of  the  full  offline  reconstruction  code.  This  means  that  for 
example  that  fully  reconstructed  3-dimensional  tracks  are  available  to  the  trigger 
decision.  The  Level-3  output  rate  is  ~  75  Hz  and  accepted  events  are  written  to  tape  in 
eight  separate  'streams',  sorted  by  the  Consumer-Server  Logger  (CSL). 


CHAPTER  3 
JETS  AT  CDF 

We  have  studied  the  "underlying  event"  in  the  Run  2  jet  trigger  data  samples  using 
the  direction  of  the  leading  calorimeter  jet  (Midpoint,  R  =  0.7,  fmerge  =  0.75)  to  isolate 
regions  of  ri-(|)  space  that  are  sensitive  to  the  "underlying  event". 

Hadronization  of  the  outgoing  partons  forms  the  jets  we  see  experimentally.  Jet 
algorithms  are  employed  to  map  data  onto  jets  with  the  idea  that  theses  jets  are  surrogates 
for  the  underlying  energetic  partons.  In  our  theoretical  picture,  the  partons  produced  by 
the  hard  scattering  process  evolve  approximately  within  a  narrow  cone  based  on  the 
parton  showering  and  hadronization  models.  We  use  the  Midpoint  [47-49]  jet  algorithm 
in  which  the  properties  of  the  constituents  of  the  jet(J)  of  cone  radius  R  are  defined  by  the 
following  equations 

keJ:{0,-(^,f+{y,-yjf<R'  3-1 

keJ     ^T,J  ksJ     '-'T.J 

^T.J   =zl^T,k  •  ^'^ 

k€j 

The  MidPoint  cone  algorithm  is  based  on  the  so-called  "Snowmass  Algorithm"  which 
defines  both  the  stability  conditions  and  the  properties  of  the  jets  [50]. 

The  basic  jet  cone  idea  is  that  the  constituents  are  nearby  each  other  in  simple 
geometric  fashion.   That  is,  the  3-momenta  of  the  hadrons  or  partons  lie  within  a  cone 

defined  by  a  circle  in  the  angular  variables  (y,(p),  where  y  =  —ln[(£' +/?,)/(£- p  J]  is 
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the  true  rapidity  and  (p  is  the  azimuthal  angle.  A  stable  jet  cone  has  the  property  that  the 
geometric  center  of  the  cone  coincides  with  the  Ej  weighted  centriod  of  the  particles  in 
the  cone.  Jet  algorithms  involve  two  distinct  steps.  The  first  is  to  identify  the 
"constituents"  that  comprise  the  stable  cone  that  is  the  jet.  The  second  involves 
constructing  the  kinematic  properties  that  characterize  the  jet.  However,  in  practice  the 
experimental  implementation  of  the  cone  algorithm  is  more  complicated. 

It  was  imagined  that  the  entire  particle/tower  list  of  each  event  would  be  searched 
for  sets  of  final  state  particle/towers  which  satisfy  the  stability  conditions.  In  practice  this 
is  not  possible  because  of  limited  computing  resources  so  a  number  of  compromises  have 
had  to  be  made. 

The  Midpoint  algorithm  starts  with  an  Ej  ordered  list  of  seed  towers  (Etseed  > 
threshold),  and  forms  "protojets"  from  every  stable  cone  iterated  around  a  seed  tower. 
Seed  Towers  are  simply  calorimeter  towers  in  which  the  energy  deposition  exceeds  a 
certain  predefined  limit  (usually  set  to  1  GeV)  which  is  larger  than  the  limit  defined  to 
include  a  tower  in  a  jet  (typically  0.3  GeV).  A  search  for  new  protojets  is  carried  out 
about  the  Midpoints  in  (y,(p)  between  all  pairs  of  protojets  with  AR  <  2xRcone-  The 
Midpoint  algorithm  includes  an  iterative  splitting/merging  process  applied  to  the  Pt 
ordered  list  of  jets  to  assign  each  particle  to  only  one  jet.  Two  jets  are  merged  if  the 
lower  Pt  jet  shares  greater  than  fmerge  of  its  total  Pt  with  the  higher  Pt  jet.  Otherwise  the 
two  jets  are  split  and  the  individual  particles  are  assigned  to  the  closest  jet  centriod. 


CHAPTER  4 
MONTE-CARLO  GENERATION  AND  CORRECTION  FACTORS 

4. 1  Monte-Carlo  Generation 

In  this  analysis  the  data  are  corrected  back  to  the  particle  level  using  PYTHIA  Tune 

A  [22,  23].   The  corrected  data  are  then  compared  with  the  particle  level  predictions  of 

PYTHIA  Tune  A  and  HERWIG  (i.e.  generator  level)  at  1.96  TeV.     PYTHIA  Tune  A 

(5.3.3nt)  was  generated  with  the  minimum  PrChard)  values  shown  in  Table  4-1  and 

HERWIG  (5.3.3nt)  was  generated  with  the  minimum  PiChard)  values  shown  in  Table  4-2. 

Stntuples  (5.3.3nt  dev242)  were  created  for  the  QCD  group  by  Anwar  Bhatti,  Ken 

Hatakeyama,  and  Craig  Group. 

Table  4-1:  PYTHIA  Tune  A  (5.3.3nt)  at  1.96  TeV. 


PxChard)  minimum 

Events 

0  GeV/c 

3,093,106 

10  GeV/c 

1,039,093 

18GeV/c 

4,285,687 

40  GeV/c 

4,228,873 

60  GeV/c 

992,087 

90  GeV/c 

1,497,108 

120GeV/c 

2,068,377 

150GeV/c 

1,488,786 

200  GeV/c 

1,042,280 

300  GeV/c 

1,045,314 

400  GeV/c 

1,043,634 

Total 

21,824,345 
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Table  4-2:  HERWIG  (5.3.3nt)  at  1.96  TeV. 


PxChard)  minimum 

Events 

3  GeV/c 

1,014,070 

lOGeV/c 

1,018,974 

18GeV/c 

5,001,261 

40  GeV/c 

5,071,205 

60  GeV/c 

1,044,202 

90  GeV/c 

2,057,661 

120GeV/c 

2,035,473 

150  GeV/c 

1,922,568 

200  GeV/c 

968,906 

300  GeV/c 

885,867 

400  GeV/c 

858,936 

Total 

21,879,123 

Smooth  curves  have  been  drawn  through  the  QCD  Monte-Carlo  predictions  to  aid 
in  comparing  the  theory  with  the  data.  Fig.  4-1  shows  an  example  of  the  fits  to  the 
Monte-Carlo  results. 

4.2  Correcting  the  Data  to  the  Particle  Level 

We  consider  two  methods  for  correcting  the  data  from  the  detector  level  to  the 
particle  level.  Method  1  is  a  "one-step"  method  in  which  PYTHIA  Tune  A  and 
HERWIG  are  used  to  calculate  the  observables  in  Table  4-3  at  the  particle  level  (in  bins 
of  particle  jet#l  Pj  "G£7V")  and  at  the  detector  level  (in  bins  of  calorimeter  jet#l  Pj 
(uncorrected)  "CDFS/Af' ).     The  detector  level  data,  in  bins  of  calorimeter  jet#l  Pj 
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(uncorrected),  are  corrected  by  multiplying  by  a  QCD  Monte-Carlo  correction  factor, 

GEN/CDFSIM,  as  described  in  Table  4-4. 

Table  4-3:  Observables  examined  in  the  "transverse"  region  (see  Fig.  6-2)  as  they  are 
defined  at  the  particle  level  and  the  detector  level.  Charged  tracks  are  considered  "good" 
if  they  pass  the  selection  criterion  given  in  Table  5-2.  The  mean  charged  particle  <pT> 
and  the  charged  fraction  PTsum/ETsum  are  constructed  on  and  event-by-event  basis  and 
then  averaged  over  the  events.  There  is  one  PTmax  per  event  with  PTmax  =  0  if  there  are 
no  charged  partic  es. 


Observable 


Particle  Level 


Detector  level 


dNchg/drid(|) 


Number  of  charged  particles 

per  unit  y\-<^ 

(PT->0.5GeV/c,  h|<l) 


Number  of  "good"  charged  tracks 

per  unit  r\-(^ 

(l>r>0.5GeV/c,hl<l) 


dPTsum/drid(t) 


Scalar  pr  sum  of  charged 

particles  per  unit  r\-(^ 
(pr>0.5GeV/c,  |ti|<1) 


Scalar  pr  sum  of  "good"  charged  tracks 

per  unit  r|-(j) 

(pr>0.5GeV/c,  |ril<l) 


<Pt> 


Average  pT  of  charged  particles 
(pT>0.5GeV/c,  h|<  1) 


Maximum  pr  charged  particle 

(pr>0.5GeV/c,  |ti|<1) 

PTmax  =  0  for  no  charged 

particle 


Average  pr  of  "good"  charged  tracks 
(pr>0.5GeV/c,  hl<l) 


PTmax 


Maximum  pr  "good"  charged  tracks 

(pT>0.5GeV/c,  |tiI<1) 
PTmax  =  0  for  no  "good"  charged  track 


dEi/dridcl) 


Scalar  Ej  sum  of  all  particles 

per  unit  r[-^ 

(aIlpT,h|<l) 


Scalar  Ex  sum  of  all  calorimeter  towers 

per  unit  T|-(|) 
(ET>0.1GeV,  hl<l) 


PTsum/ETsum 


Scalar  pr  sum  of  charged 

particles 

(pr>0.5GeV/c,  |ti|<1) 

divided  by  the  scalar  Ex  sum  of 

all  particles  (all  px,  |t||  <  1) 


Scalar  px  sum  of  "good"  charged  tracks 

(pT>0.5GeV/c,  |tiI<1) 

divided  by  the  scalar  Ex  sum  of 

calorimeter  towers 

(Ex > 0.1  GeV,  fTi|<l) 


Method  2  is  a  "two-step"  method.  First  PYTHIA  Tune  A  is  used  to  correct  the  Pj 
of  the  leading  calorimeter  jet.  This  is  done  by  comparing  the  matching  leading  particle 
jet  with  the  leading  calorimeter  jet.  Then  PYTHIA  Tune  A  is  used  to  calculate  the 
observables  in  Table  4-3  at  the  particle  level  in  bins  of  particle  jet#l  Px  (GEN)  and  at  the 
detector  level  in  bins  of  calorimeter  jet#l  Px  (corrected)  (CDFSIMcor).  The  detector 
level  data  in  bins  of  calorimeter  jet#l  Px  (corrected)  are  corrected  by  multiplying  by  the 
QCD  Monte-Carlo  correction  factor,  GEN/CDFSIMcor.     If  the  QCD  Monte-Carlo 
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described  the  data  perfectly  and  the  detector  simulation  was  exact  then  method  1  and 

method  2  would  yield  the  same  result.  Differences  between  the  two  methods  can  be  used 

as  a  measure  of  the  systematic  uncertainty  in  correcting  the  data  to  the  particle  level. 

Table  4-4:  Correction  factors  for  Method  1.  PYTHIA  Tune  A  and  HERWIG  are  used  to 
calculate  the  observables  in  Table  4-3  at  the  particle  level  in  bins  of  particle  jet#l  Pj 
(GEN)  and  at  the  detector  level  in  bins  of  calorimeter  jet#l  Pt  (uncorrected).  The 
detector  level  data  in  bins  of  calorimeter  jet#l  PT  (uncorrected)  are  corrected  by 
multiplying  by  QCD  Monte-Cario  factor,  GEN/CDFSIM. 


Particle  Level 
Observable 

Detector  Level 
Observable 

Response  Factor 

Correction 
Factor 

GEN  =  Particle 
Jet#l  Pt  Bin 

CDFSIM  =  Calorimeter 
Jet#l  Pt  Bin  (uncorrected) 

CDFSIM/GEN 

GEN/CDFSIM 

Fig.  4-2  and  Fig.  4-3  show  the  particle  level  predictions  from  PYTHIA  Tune  A  and 
HERWIG  for  average  density  of  particles,  dNaii/drjdct),  for  all  particles  with  |t||  <  1  in  the 
"transverse"  region  as  a  function  of  the  leading  particle  jet  Pj  for  "leading  jet"  and 
"Back-to-back"  events,  respectively  (see  figure  6-3).  It  is  interesting  to  note  that 
HERWIG  produces  more  particles  in  the  "transverse"  region  than  PYTHIA  Tune  A.  Fig. 
4-2  and  Fig.  4-3  also  shows  the  average  charged  particle  PTsum  density,  dPTsum/dridcl),  and 
the  average  charged  particle  <Pt>  for  particles  with  |t||  <  1  in  the  "transverse"  region  for 
"leading  jet"  events  as  a  function  of  the  leading  particle  jet  Pj.  It  is  clear  from  these 
comparisons  that  HERWIG  produces  more  "soft"  particles  than  PYTHIA  Tune  A  which 
will  result  in  different  "response"  factors  (see  Table  4-4)  at  low  leading  jet  Pj. 
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Figure  4-1:  Example  of  fits  to  the  QCD  Monte-Carlo  results.  Shows  the  particle  level 
predictions  at  1 .96  TeV  for  the  density  of  charged  particles,  dNchg/dr|d(l),  with  px  >  0.5 
GeV/c  and  |ti|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "leading  jet"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  particle  jet  Pt  for  PYTHIA  Tune  A  {top) 
and  HERWIG  (bottom). 
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Figure  4-2:  Particle  level  predictions  from  PYTHIA  Tune  A  and  HER  WIG  for  average 
density  of  particles  dNall/dr|d(l)  (top),  the  average  charged  particle  PTsum  density, 
dPTsum/dr|d(t)  (middle),  and  the  average  charged  particle  <Pt>  (bottom)  for  particles  with 
|t||  <  1  in  the  "transverse"  region  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function 
of  the  leading  particle  jet  Pt. 


44 


Transverse"  Particle  Density:  dN/dridtt) 


5.0 


5  4.0 


.u  3.0 


%   2.0 


1.0 


0.0 


CDF  Run  2  Preliminary 

generator  lever  theory 

"Back-to-Back* 

*    _    _     _    _ 

HWj 



/^ 

PYTuneA  1 

MIdPoint  R  =  0.  /  |Ti(jet)|  <  2 

1.96  TeV 

i 1 1 ^ 

All  Particles  (hi |<1.0,  all  PT) 

M 1 1 1 1 

50    100    150    200    250    300    350 

PT(particle]et#1)  (GeV/c) 


400        450 


500 


1.8 


'Transverse"  PTsum  Density:  dPT/dTid(|) 


> 

o 

a 
« 

c 
o 
a 

E 


1.2 


•0.6 


0.0 


CDF  Run  2  Preliminary 

generator  level  theory 

1.96  TeV 


"Back-to-Back" 


Midpoint  R  =  0.7  |T|Oet)|  <  2 
Charged  Particles  (|i)|<1.0,  all  PT) 


-h 


-+- 


-+- 


-+- 


50  100         150        200        250        300        350 

PT(partlclejet#1)    (GeV/c) 


400        450 


500 


'Transverse"  Average  Ctiarged  PT 


1.25 


>  1.00 


0.75 


0.50 


0.25 


CDF  Run  2  Preliminary 

generator  level  theory 

f.96TeV 


Charged  Particles  (|ii|<1.0,  all  PT) 
MIdPoint  R  =  0.7  |i)Uet)|  <  2 

"Back-to-Back" 


PY  Tune  A 


[hw] 


—I 1 1 1 1 1 1 1 1 

50         100        150        200        250        300        350        400        450        500 
PT(partlclejet#1)  (GeV/c) 


Figure  4-3:  Particle  level  predictions  from  PYTHIA  Tune  A  and  HER  WIG  for  average 
density  of  particles  dNall/dridcl)  (top),  the  average  charged  particle  PTsum  density 
dPTsum/dr|d(t)  (middle),  and  the  average  charged  particle  <pT>  (bottom)  for  particles  with 
|ti|  <  1  in  the  "transverse"  region  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a 
function  of  the  leading  jet  Pj. 
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Figs.  4-4  thru  4-22  show  the  "response"  factors  (see  Table  4-4)  from  PYTHIA  Tune 
A  and  HERWIG  for  the  observables  in  Table  4-3  for  "leading  jet"  and  "back-to-back" 
events  as  a  function  of  the  leading  jet  Px.  HERWIG  and  PYTHIA  Tune  A  produce 
similar  "response"  factors  for  leading  jet  Pt  greater  than  50  GeV/c,  but  for  lower  leading 
jet  Pt  they  are  quite  different.  This  will  result  in  large  systematic  errors  on  the  corrected 
observables  in  Table  4-3  at  low  leading  jet  Pt. 

Fig.  4-23  shows  the  leading  jet  Pt  correction  used  in  method  2  for  "leading  jet" 
events.  Figs.  4-24  shows  the  method  2  "response"  factors  from  PYTHIA  Tune  A  for 
some  of  the  observables  in  Table  4-3  for  "leading  jet"  events  as  a  function  of  the  leading 
jet  Pt.  The  observable  in  Table  4-3  do  not  depend  strongly  on  the  leading  jet  Pt  and 
hence  the  method  1  and  method  2  correction  factors  are  similar.  This  can  be  seen  in  Fig. 
4-25  which  compares  the  method  1  response  factors  versus  the  leading  jet  Pt 
(uncorrected)  with  the  method  2  response  factors  versus  the  leading  jet  Pt  (corrected) 
from  PYTHIA  Tune  A.  The  method  2  correction  factors  (1 /response  factor)  are  apphed 
data  after  correcting  the  leading  jet  Pt,  while  the  method  1  correction  factors  are  applied 
to  the  data  without  correcting  the  leading  jet  Pt. 

Method  1  can  be  easily  applied  to  both  the  "leading  jet"  and  "back-to-back"  events. 
In  "back-to-back"  events,  method  1  corrects  for  calorimeter  response  for  jet#l,  jet#2,  and 
jet#3  in  one  step.  Figures  4-4  through  4-22  show  that  the  response  factors  for  "back-to- 
back"  events  are  different  from  those  of  the  "leading  jet"  events.  The  primary  source  of 
the  difference  is  due  to  the  requirement  that  PT(jet#3)  <  15  GeV/c  for  "back-to-back" 
events  and  the  "back-to-back"  correction  factors  are  correcting  for  the  calorimeter 
response  for  jet#3.  In  order  to  apply  method  2  to  the  "back-to-back"  events  we  would 
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Figure  4-4:  Method  1  response  factors  for  the  density  of  charged  particles,  dNchg/dTidcj), 
with  pT  >  0.5  GeV/c  and  |ti|  <  1  in  the  "transMAX"  region  for  "leading  jet"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt.  Shows  the  particle  level  prediction 
(GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result  (CDFSEM)  versus 
the  leading  calorimeter  jet  Pj  (uncorrected)  with  |r|(jet#l)|  <  2  for  PYTHIA  Tune  A  (top) 
and  HER  WIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pj  (i.e.  response  factor). 
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have  to  first  correct  the  Pt  of  jet#l,  jet#2,  and  jet#3.  Since  the  observables  in  Table  4-3 
do  not  depend  strongly  on  the  Pt  of  jet#l,  jet#2,  and  jet#3,  it  is  much  easier  to  use 
method  1  for  both  "leading  jet"  and  "back-to-back"  events.  We  will  use  the  differences 
between  method  1  and  method  2  in  "leading  jet"  events  as  a  measure  of  the  systematic 
uncertainty  in  correcting  to  the  particle  level. 
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Figure  4-5:  Method  1  response  factors  for  the  density  of  charged  particles,  dNchg/drid(j), 
with  Pt  >  0.5  GeV/c  and  |ti|  <  1  in  the  "transMIN"  region  for  "leading  jet"  events  defined 
in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt.  Shows  the  particle  level  prediction  (GEN) 
versus  the  leading  particle  jet  Pj  and  the  detector  level  result  (CDFSIM)  versus  the 
leading  calorimeter  jet  Pt  (uncorrected)  with  |ri(jet#l)|  <  2  for  PYTHIA  Tune  A  (top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-6:  Method  1  response  factors  for  the  PTsum  density  of  charged  particles, 
dPTsum/dridct),  with  pj  >  0.5  GeV/c  and  |ri|  <  1  in  the  "transMAX"  region  for  "leading 
jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Ft.  Shows  the  particle  level 
prediction  (GEN)  versus  the  leading  particle  jet  Ft  and  the  detector  level  result 
(CDFSIM)  versus  the  leading  calorimeter  jet  Ft  (uncorrected)  with  |T|(jet#l)|  <  2  for 
FYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level 
to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Fj  (i.e.  response  factor). 
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Figure  4-7:  Method  1  response  factors  for  the  PTsum  density  of  charged  particles, 
dPTsum/dTid(t),  with  pj  >  0.5  GeV/c  and  |ri|  <  1  in  the  "transMESf"  region  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt-  Shows  the  particle  level 
prediction  (GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result 
(CDFSIM)  versus  the  leading  calorimeter  jet  Pt  (uncorrected)  with  |r|(jet#l)|  <  2  for 
PYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level 
to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-8:  Method  1  response  factors  for  the  average  <Pt>  of  charged  particles  with  px  > 
0.5  GeV/c  and  |t||  <  1  in  the  "transverse"  region  for  "leading  jet"  events  defined  in  Fig.  6- 
3  as  a  function  of  the  leading  jet  Px.  Shows  the  particle  level  prediction  (GEN)  versus  the 
leading  particle  jet  Px  and  the  detector  level  result  (CDFSIM)  versus  the  leading 
calorimeter  jet  Px  (uncorrected)  with  |Ti(jet#l)|  <  2  for  PYTHIA  Tune  A  (top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Px  (i.e.  response  factor). 
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Figure  4-9:  Method  1  response  factors  for  the  average  maximum  px,  PTmax,  for  charged 
particles  with  px  >  0.5  GeV/c  and  |r||  <  1  in  the  "transverse"  region  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px.  Shows  the  particle  level 
prediction  (GEN)  versus  the  leading  particle  jet  Px  and  the  detector  level  result 
(CDFSIM)  versus  the  leading  calorimeter  jet  Px  (uncorrected)  with  |ri(jet#l)|  <  2  for 
PYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level 
to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Px  (i.e.  response  factor). 


53 


5.0 


■TransMAX"  ETsum  Density:  dET/dTid^) 


0.0 


CDF  Run  2  Preliminary 

PYTHIATuneA 


beading  def^ 

Midpoint  R  =  0.7  |i|(iet)|  <  2 


Towers  (|t)|<1.0,  ET>0.1  GeV) 


50         100        150        200        250        300        350        400        450 
PT(jet#1  uncorrected)  or  PT(particle  jet#1)  (GeV/c) 


500 


TransMAX"  ETsum  Density:  dET/diidi|) 


5.0 


,4.0 


2.0 


£  1-0 


0.0 


CDF  Run  2  Preliminary 

HERWIG 

^....■^—^^■^^^ 

'"^                          "^   GEN    1 

1.96  TeV^'''^'^ 

'"' "^                                                    ~~ 

COFSIM     1 

•Leading  Jet" 
Midpoint  R  =  0.7  h|(iet)|  <  2 

/        y^      1  CDFSiMAjEN  | 

1 \ \ 

Charged  Particles  (|ti|<1  .0,  PT>0.5  GeV/c) 

\ \ \ \ 1 1 

50         100        150        200        250        300        350        400        450 
PT(jet#1  uncorrected)  or  PT(particlejet#1)  (GeV/c) 


500 


TransMAX"  ETsum  Density:  CDFSIM/GEn] 


1.0 

0.9 

I  0.8 

s 

in 


iO.6 

cc 


0.5 
0.4 


CDF  Run  2  Preliminary                                           -lyeadingJei" 

1.96  TeV             ^j;-i>^ 

.„^^ 

i/^^^"^ ^ ^*-— ^ 

K ^^— -— 

1                  1                                                                            MiHDnin*  R        H  7  Ut/iat\I  ^  9 

Ciiarged  Particles  (|iil<1.0,  PT>0.5  GeV/c) 

1 1 1 \ \ 1 1 1 1 

50         100        150        200        250        300        350 
PT(jet#1  uncorrected)  (GeV/c) 


400        450 


500 


Figure  4-10:  Method  1  response  factors  for  the  ETsum  density  of  all  particles,  dEj/dridct), 
with  |ri|  <  1  in  the  "transMAX"  regions  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a 
function  of  the  leading  jet  Pj.  Shows  the  particle  level  prediction  (GEN)  versus  the 
leading  particle  jet  Pt  and  the  detector  level  result  (CDFSIM)  versus  the  leading 
calorimeter  jet  Pt  (uncorrected)  with  |T|(jet#l)|  <  2  for  PYTHIA  Tune  A  {top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-11:  Method  1  response  factors  for  the  ETsum  density  of  all  particles,  dEj/dridct), 
with  |r||  <  1  in  the  "transMIN"  regions  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a 
function  of  the  leading  jet  Pt.  Shows  the  particle  level  prediction  (GEN)  versus  the 
leading  particle  jet  Pj  and  the  detector  level  result  (CDFSIM)  versus  the  leading 
calorimeter  jet  Pt  (uncorrected)  with  |Ti(iet#l)|  <  2  for  PYTHIA  Tune  A  (top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pj(i.e.  response  factor). 
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Figure  4-12:  Method  1  response  factors  for  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Pj,  where  PTsum  includes  charged  particles  with  px  >  0.5  GeV/c  and  |t||  <  1 
and  the  ETsum  includes  all  particles  with  |ri|  <  1.  Shows  the  particle  level  prediction 
(GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result  (CDFSIM)  versus 
the  leading  calorimeter  jet  Pt  (uncorrected)  with  |Ti(jet#l)|  <  2  for  PYTHIA  Tune  A  (top) 
and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  PT(i.e.  response  factor). 


56 


1.2 -| 

1.1 

z 
111 

s  1.0 

<n 
u. 

o 

« 
0.8 

0.7 

t 

"Transverse"  Charged  Particle  Density:  CDFSIM/GEn| 

1         "Leading  Jet"                                        CDF  Run  2  Preliminary 

B                                                                                                                PYTHIATuneA 

1                                 |Tran.MAx|                                      ^'^'^^^ 

1 ""^^^ i  TransAVE  | 

r ^^            t-^^ 

.J^_^_^____^^_^^^^^^^^                                 _^ 

1   TransMIN   1 

i                   r 
MidPoint  R  =  0.7  InOet)!  <  2                    Charged  Particles  (|ii|<1 .0,  PT>0.5  GeV/c) 

)           50         100        150        200        250        300        350        400        450        500 
PT(jet#1  uncorrected)  (GeV/c) 

'Transverse"  Charged  PTsum  Density:  CDFSIIWGENJ 


1.2 

1.1 


O  0.9 
u 


0.8 
0.7 


1                                                                    CDF  Run  2  Preliminary 

1                                                                                                                PYTHIA  Tune  A 

1 

1.96  TeV 

'l     ^-^ 

I  TransMAX 

fc; 

\  \ 

■««v,^ 

^^ 

y                     1  TransAVE  | 

:S=a^^^^^ 

tMIN   1 





1  Tram 

MIdPotnt  R 

1 

=  0.7|iiaet)|< 

— 1 V- 

2 
1 

Charged  Particles  (|t)|<1.0,  PT>O.S  GeV/c) 

H 1 1 1 1 

50         100        150        200        250        300        350        400 
PT(jet#1  uncorrected)  (GeV/c) 


450        500 


"Transverse"  ETsum  Density:  CDFSII«/GEn| 


1.0 


0.9- 


COF  Run  2  Preliminary 

PYTHIATuneA 


"Leading  Jet" 


1.96  TeV 


50         100        150        200        250        300        350        400        450        500 
PT(iet#1  uncorrected)  (GeV/c) 


Figure  4-13:  Shows  the  ratio  of  the  detector  level  to  the  particle  level,  CDFSIM/GEN, 
versus  the  leading  jet  Pj  (method  1  response  factors)  for  PYTHIA  Tune  A  for  the 
"transMAX",  "transMIN",  and  "transverse"  regions  for  "leading  jet"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Pt.  Shows  the  density  of  charged  particles 
dNchg/dr|d(t)  with  pi  >  0.5  GeV/c  and  hi  <  1  {top),  the  PTsum  density  of  charged 
particles  dPTsum/drid(l)  with  px  >  0.5  QtSIc  and  |r||  <  1  (middle),  and  ETsum  density  of 
all  particles  dEj/dr\d(^  with  hi  <  1  (bottom). 
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Figure  4-14:  Method  1  response  factors  for  the  density  of  charged  particles,  dNchg/dridcj), 
with  pt  >  0.5  GeV/c  and  |r||  <  1  in  the  "transMAX"  region  for  "back-to-back"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt-  Shows  the  particle  level  prediction 
(GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result  (CDFSIM)  versus 
the  leading  calorimeter  jet  Pt  (uncorrected)  with  |r|(jet#l)|  <  2  for  PYTHIA  Tune  A  (top) 
and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-15:  Method  1  response  factors  for  the  density  of  charged  particles,  dNchg/dT|d(|), 
with  pt  >  0.5  GeV/c  and  |r||  <  1  in  the  "transMIN"  region  for  "back-to-back"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj.  Shows  the  particle  level  prediction 
(GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result  (CDFSIM)  versus 
the  leading  calorimeter  jet  Pj  (uncorrected)  with  ItiOetttl)]  <  2  for  PYTHIA  Tune  A  {top) 
and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-16:  Method  1  response  factors  for  the  PTsum  density  of  charged  particles, 
dPTsum/dridcl),  with  px  >  0.5  GeV/c  and  |ti|  <  1  in  the  "transMAX"  region  for  "back-to- 
back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt.  Shows  the  particle 
level  prediction  (GEN)  versus  the  leading  particle  jet  Pj  and  the  detector  level  result 
(CDFSIM)  versus  the  leading  calorimeter  jet  Pj  (uncorrected)  with  |Ti(jet#l)|  <  2  for 
PYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level 
to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Py  (i.e.  response  factor). 
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Figure  4-17:  Method  1  response  factors  for  the  PTsum  density  of  charged  particles, 
dPTsum/dr|d(t),  with  px  >  0.5  GeV/c  and  |r||  <  1  in  the  "transMIN"  region  for  "back-to- 
back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj.  Shows  the  particle 
level  prediction  (GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result 
(CDFSIM)  versus  the  leading  calorimeter  jet  Pt  (uncorrected)  with  |Ti(jet#l)|  <  2  for 
PYTHIA  Tune  A  (top)  and  HERWIG  {middle).  Also  shows  the  ratio  of  the  detector  level 
to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Pj  [i.e.  response  factor). 
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Figure  4-18:  Method  1  response  factors  for  the  average  <Pt>  of  charged  particles  with  pr 
>  0.5  GeV/c  and  |ri|  <  1  in  the  "transverse"  region  for  "back-to-back"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Py.  Shows  the  particle  level  prediction  (GEN) 
versus  the  leading  particle  jet  Pt  and  the  detector  level  result  (CDFSIM)  versus  the 
leading  calorimeter  jet  Pt  (uncorrected)  with  |r|(jet#l)|  <  2  for  PYTHIA  Tune  A  (top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-19:  Method  1  response  factors  for  the  average  maximum  pr,  PTmax,  for  charged 
particles  with  pi  >  0.5  GeV/c  and  ItjI  <  1  in  the  "transverse"  region  for  "back-to-back" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj.  Shows  the  particle  level 
prediction  (GEN)  versus  the  leading  particle  jet  Pt  and  the  detector  level  result 
(CDFSIM)  versus  the  leading  calorimeter  jet  Pt  (uncorrected)  with  |ri(jet#l)|  <  2  for 
PYTHIA  Tune  A  (top)  and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level 
to  the  particle  level,  CDFSIM/GEN,  versus  the  leading  jet  Pt  (i.e.  response  factor). 
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Figure  4-20:  Method  1  response  factors  for  the  ETsum  density  of  all  particles,  dEx/dridcj), 
with  |t||  <  1  in  the  "transMAX"  regions  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a 
function  of  the  leading  jet  Pj.  Shows  the  particle  level  prediction  (GEN)  versus  the 
leading  particle  jet  Pj  and  the  detector  level  result  (CDFSIM)  versus  the  leading 
calorimeter  jet  Pt  (uncorrected)  with  |r|(iet#l)|  <  2  for  PYTHIA  Tune  A  {top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Pj  (i.e.  response  factor). 
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Figure  4-21:  Method  1  response  factors  for  the  ETsum  density  of  all  particles,  dEt/dridcl), 
with  |ti|  <  1  in  the  "transMIN"  regions  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a 
function  of  the  leading  jet  Pj.  Shows  the  particle  level  prediction  (GEN)  versus  the 
leading  particle  jet  Pt  and  the  detector  level  result  (CDFSIM)  versus  the  leading 
calorimeter  jet  Pt  (uncorrected)  with  |Ti(iet#l)|  <  2  for  PYTHIA  Tune  A  (top)  and 
HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  PT(i.e.  response  factor). 
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Figure  4-22:  Method  1  response  factors  for  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Pt,  where  PTsum  includes  charged  particles  with  px  >  0.5  GeV/c  and  hi  <  1 
and  the  ETsum  includes  all  particles  with  |ti|  <  1.  Shows  the  particle  level  prediction 
(GEN)  versus  the  leading  particle  jet  Py  and  the  detector  level  result  (CDFSIM)  versus 
the  leading  calorimeter  jet  Fj  (uncorrected)  with  |r|(jet#l)|  <  2  for  PYTHIA  Tune  A  (top) 
and  HERWIG  (middle).  Also  shows  the  ratio  of  the  detector  level  to  the  particle  level, 
CDFSIM/GEN,  versus  the  leading  jet  Px  (i.e.  response  factor). 
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Figure  4-23:  Leading  jet  Px  correction  used  in  method  2  for  "leading  jet"  events.  Shows 
the  difference  in  the  observed  leading  jet  Pt  at  the  detector  level  (JL.e.  in  the  calorimeter) 
compared  with  the  true  Pj  (j.e.  corrected)  of  matched  leading  particle  jets  using  PYTHIA 
Tune  A  and  HERWIG. 
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Figure  4-24:  Method  2  response  factors  from  PYTHIA  Tune  A  for  "leading  jet"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt.  Shows  the  particle  level  prediction 
(GEN)  versus  the  leading  particle  jet  Pj  and  the  detector  level  result  (CDFSIMcor)  versus 
the  leading  calorimeter  jet  Pj  (corrected)  with  |ii(jet#l)|  <  2.  Also  shows  the  ratio  of  the 
detector  level  to  the  particle  level,  CDFSIMcor/GEN,  versus  the  leading  jet  Pt  {i.e. 
response  factor). 
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Figure  4-25:  Compares  the  method  1  response  factors  versus  the  leading  jet  Pj 
(uncorrected)  with  the  method  2  response  factors  versus  the  leading  jet  Pt  (corrected) 
from  PYTHIA  Tune  A. 


CHAPTER  5 
DATA  SELECTION  AND  SYSTEMATIC  ERRORS 

5.1  Data  Selection 
The  data  used  in  this  analysis  arise  from  the  set  of  Stntuples  created  for  the  QCD 

group  by  Anwar  Bhatti,  Ken  Hatakeyama,  and  Craig  Group  (see  Table  5-1).  Events  are 

required  to  be  on  the  "goodrun"  list  (version  7).  They  are  also  required  to  have  a  niissing 

Et  significance  less  than  5  GeV"^  and  to  have  a  sumET  <  1.5  TeV.  Except  for  the  Min- 

Bias  data  we  require  events  to  have  one  and  only  one  quality  12  vertex  with  |z|  <  60  cm. 

For  the  Min-Bias  data  we  allow  zero  or  one  quality  12  vertices.    This  only  affects  the 

observables  in  Table  4-3  for  leading  jet  Pt  below  10  GeV/c.  Above  PjOetttl)  =  10  GeV/c 

the  fraction  of  events  with  no  quality  12  vertex  is  negligible. 

Table  5-1:  Data  sets  (5.3.3nt)  and  event  selection  criterion  used  in  this  analysis 
(L  ~  380  pb"'). 


Event  Selection 

Min-Bias 

JET20 

JET50 

JET70 

JET  100 

Total  Events 

20,586,733 

30,470,383 

9,908,366 

4,641,247 

5,366,515 

"Good"  Events  (version  7) 

18,180,015 

19,835,681 

6,868,114 

3,432,992 

4,031,324 

MetSig<5GeV"^, 
sumET<1.5TeV 

18,179,280 

19,818,879 

6,785,357 

3,316,514 

3,602,989 

1  Q12ZVtx,  z  <60cm 

15,416,180 

10,851,963 

3,745,616 

1,794,739 

1,939,382 

"Leading  Jet"  ri(jet#l)  <  2 

3,712,407 

7,679,594 

3,200.065 

1,648,764 

1,884,353 

"Back-to-Back" 
PT(jet#3)  <  15  GeV/c 

2,474 

1,462,547 

878,014 

491,930 

602,256 

"Back-to-Back"/"Leading 
Jet" 

0.07% 

19.04% 

27.44% 

29.84% 

31.96% 

As  in  a  Run  1  analysis  [13]  only  charged  particles  in  the  region  px  >  0.5  GeV/c  and  |ti|  <  1 
where  the  COT  efficiency  is  high  are  considered.  Our  track  selection  criterion  shown  in 
Table  5-2  is  the  same  as  the  Run  1  analysis. 
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Table  5-2:  Track  Selection  criterion. 


Track  Selection 


COT  measured  tracks 


|z-zo|  <  2  cm 


|do|  <  1  cm 


pr>0.5GeV/c,  |ti|<1 


In  forming  the  observables  in  Table  4-3  the  five  trigger  sets  shown  in  Table  5-1  are 
pieced  together  as  shown  in  Table  5-3.   The  "looser"  trigger  set  is  used  until  it  overlaps 
the  next  trigger  set  and  then  that  trigger  set  is  used  until  it  overlaps  the  next  trigger  set  etc. 
Table  5-3:  Range  of  PT(jet#l)  used  for  each  data  set. 


Trigger  Set 


Min-Bias 


JET20 


JET50 


JET70 


JET  100 


Calorimeter  Jets 


PT(iet#l)<30GeV 


30<PTGet#l)<70GeV 


70<PT(jet#l)<95GeV 


95<PTaet#l)<130GeV 


PT(iet#l)>130GeV 


5.2  Systematic  Uncertainty 

The  systematic  uncertainty  in  correcting  to  the  particle  level  is  estimated  by 

combining  the  two  factors  shown  in  Table  5-4.  The  first  factor,  C\,  comes  from  correcting 
the  observables  in  Table  4-3  to  the  particle  level  using  method  1  and  examining  the  bin- 
by-bin  difference  between  PYTHIA  Tune  A  and  HERWIG  for  each  observable.  The 
second  factor,  02,  is  set  large  enough  to  include  the  differences  between  method  1  and 
method  2  and  pile-up  (only  affects  the  transverse  energy). 

Figure  5-1  shows  the  data  at  1.96  TeV  corrected  to  the  particle  level  using  method  1 
and  method  2.  The  open  red  squares  are  the  data  corrected  to  the  particle  level  using 
method  1  with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty 
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(see  Table  5-4).  The  black  dots  are  the  data  corrected  to  the  particle  level  using  method  2 

(with  no  errors).  The  method  2  points  lie  within  the  errors  of  the  method  1  data  points. 

Table  5-4:  The  errors  on  the  corrected  observables  in  Table  4-3  include  both  the  statistical 
error  and  the  systematic  uncertainty  (added  in  quadrature).  The  systematic  uncertainty 
consists  of  ai  and  Cj  (added  in  quadrature). 


Uncertainty 

Origin 

Oi 

Bin  by  bin  difference  between  the  data 

corrected  by  PYTHIA  Tune  A  and 

HER  WIG  using  method  1. 

02 

Difference  between  method  1  and 

method  2 

and  pile-up  and  miscellaneous 

(3%  for  charged  particle,  5%  for  energy) 
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Figure  5-1:  Data  at  1.96  TeV  corrected  to  the  particle  level  using  method  1  and  method  2 
compared  with  PYTHIA  Tune  A  and  HERWIG  at  the  particle  level.    Shows  the  density 
of  charged  particles,  dNchg/dridcj)  (top),  the  PTsum  density  of  charged  particles, 
dPTsum/drid(|)  (middle),  (px  >  0.5  GeV/c  and  hi  <  1),  and  the  ETsum  density,  dEj/dTidct) 
(bottom),  for  particles  with  |ti|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet 
Pt. 


CHAPTER  6 
DISCUSSION  OF  RESULTS 

We  study  the  behavior  of  the  charged  particles  (Pt>  0.5  GeV/c,  |r||<l)  in  the 
"underlying  event"  for  hard  scattering  pp  collisions  at  the  Tevatron  {^Js  =1.96  TeV). 
These  results  are  compared  to  QCD  Monte-Carlo  models  (PYTHIA  Tune  A  and 
HER  WIG)  that  simulate  pp  collisions.  The  topology  of  a  hard  scattering  jet  event  serves 
to  define  a  frame  of  reference.  The  direction  of  the  leading  calorimeter  jet,  jet#l,  is  used 
to  define  correlations  in  the  azimuthal  angle,  Acj).  The  angle  A(|)  =  (])  -  ^jem\  is  the  relative 
angle  between  a  charged  particle  and  the  leading  jet  direction.  Figure  6-1  shows  how  we 
partition  J]-^  space  on  an  event-by-event  basis. 

In  figure  6-2  we  define  "transverse  1"  and  "transverse  2"  so  that  we  may  compare 
these  regions  on  an  event-by-event  basis.  This  allows  us  the  flexibility  to  redefine  how 
the  regions  are  characterized  for  different  analyses.  Here  we  will  refer  to  events  in  which 
there  are  no  restrictions  placed  on  the  second  highest  Ej  jet,  jet#2,  as  "leading  jet"  events. 
Additionally,  we  define  a  subset  of  these  as  "back-to-back"  events  in  which  the  leading 
two  jets  are  nearly  "back-to-back"  (A(t)i2>150°)  and  with  Ej  (jet#2)/ETGet#l)  >  0.8. 
Within  this  subset,  the  hard  component  of  the  "underlying  event"  should  be  suppressed, 
thus  increasing  the  sensitivity  of  the  "transverse"  region  to  the  "beam-beam  remnant"  and 
multiple  parton  scattering  components. 
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Jet#l 
Direction 


Figure  6-1:  Illustration  of  correlations  in  azimuthal  angle  A(|)  relative  to  the  direction  of 
the  leading  jet  (MidPoint,  R  =  0.7,  fmerge  =  0.75)  in  the  event,  jet#l.  The  angle  A(t)  =  (])  - 
(jyetttl  is  the  relative  azimuthal  angle  between  charged  particles  (or  calorimeter  towers) 
and  the  direction  of  jet#l.  The  "toward"  region  is  defined  by  \A(^\  <  60°  and  hi  <  1,  while 
the  "away"  region  is  |A<|)|  >  120°  and  |r||  <  1.  The  "transverse"  region  is  defined  by  60°  < 
|A(t)|  <  120°  and  hi  <  1.  Each  of  the  three  regions  "toward",  "transverse",  and  "away"  and 
has  an  overall  area  in  r[-^  space  of  AriA(})  =  4n/3.  We  examine  charged  particles  in  the 
range  px  >  0.5  GeV/c  and  hj  <  1  and  calorimeter  towers  with  Et  >  0.1  GeV  and  h|  <  1, 
but  allow  the  leading  jet  to  be  in  the  region  h(jet#l)|  <  2. 

As  in  figure  6-3,  we  use  the  direction  of  jet#l  to  define  the  two  "transverse" 
regions,  60°  <  |A(j)|  <  120°.  On  an  event-by-event  basis,  we  define  the  "transMAX" 
("transMIN")  to  be  the  maximum  (minimum)  of  these  two  regions.  When  looking  at 
muldplicities  MAX  and  MIN  refer  to  the  number  of  charged  particles,  but  when 
considering  Pjsum  then  MAX  and  MIN  refer  to  the  scalar  Pxsum  of  the  charged  particles. 

Figure  6-4  illustrates  the  topology  of  a  pp  collision  in  which  a  "hard"  parton- 
parton  scattering  has  occurred.  The  contribution  from  the  hardest  initial  or  final-state 
radiation  should  be  found  in  the  "transMAX"  region.  Since  both  regions  should  receive 
"beam-beam  remnant"  contributions,  the  difference  between  "transMAX"  and 
"transMIN"  should  be  very  sensitive  to  the  "hard  scattering"  component  of  the 
"underlying  event". 


75 


Jet  #1  Direction 

A(|> 


Jet  #1  Direction 
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Figure  6-2:  Dlustration  of  correlations  in  azimuthal  angle  Acj)  relative  to  the  direction  of 
the  leading  jet  (highest  Pt  jet)  in  the  event,  jet#l .  The  angle  Acj)  =  (|)  -  (t)jet#l  is  the  relative 
azimuthal  angle  between  charged  particles  and  the  direction  of  jet#l.  The  "toward" 
region  is  defined  by  |A(t)|  <  60°  and  |ri|  <  1,  while  the  "away"  region  is  \A^\  >  120°  and  hi 
<  1 .  The  two  "transverse"  regions  60°  <  Acjx  1 20°  and  60°  <  -A^  <  120°  are  referred  to  as 
"transverse  1"  and  "transverse  2".  Each  of  the  two  "transverse"  regions  have  an  area  in 
r[-(^  space  of  AtiA(1)  =  47t/6.  The  overall  "transverse"  region  defined  in  Fig.  3  corresponds 
to  combining  the  "transverse  1"  and  "transverse  2"  regions.  Events  in  which  there  are  no 
restrictions  placed  on  the  on  the  second  highest  ?j  jet,  jet#2,  are  referred  to  as  "leading 
jet"  events  (left).  Events  with  at  least  two  jets  where  the  leading  two  jets  are  nearly 
"back-to-back"  (Ai^n  >  150°)  with  PT0et#2)/PT0et#l)  >  0.8  and  PjQem)  <  15  GeV/c  are 
referred  to  as  "back-to-back"  events  {right). 
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Jet  #2  Direction 


Fig.  6-3:  Illustration  of  correlations  in  azimuthal  angle  Acj)  relative  to  the  direction  of  the 
leading  jet  (highest  Pt  jet)  in  the  event,  jet#l  for  "leading  jet"  events  (left)  and  "back-to- 
back"  events  (right)  as  defined  in  Fig.  6-2.  The  angle  Acj)  =  (])  -  (t)jet#l  is  the  relative 
azimuthal  angle  between  charged  particles  (or  calorimeter  towers)  and  the  direction  of 
jet#l.  On  an  event  by  event  basis,  we  define  "transMAX"  ("transMIN")  to  be  the 
maximum  (minimum)  of  the  two  "transverse"  regions,  60°  <  Acjx  120°  and  60°  <  -A(^  < 
120°.  "TransMAX"  and  "transMIN"  each  have  an  area  in  V[-(^  space  of  ArjAcj)  =  47t/6. 
The  overall  "transverse"  region  defined  in  Fig.  3  includes  both  the  "transMAX"  and  the 
"transMIN"  region. 
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Figure  6-4:  Illustration  of  the  topology  of  a  proton-antiproton  collision  in  which  a  "hard" 
parton-parton  collision  has  occurred.  The  "toward"  region  as  defined  in  Fig.  6-1  contains 
the  leading  "jet",  while  the  "away"  region,  on  the  average,  contains  the  "away-side" 
"jet".  The  "transverse"  region  is  perpendicular  to  the  plane  of  the  hard  2-to-2  scattering 
and  is  very  sensitive  to  the  "underlying  event".  For  events  with  large  initial  or  final-state 
radiation  the  "transMAX"  region  defined  in  Fig.6-3  would  contain  the  third  jet  while  both 
the  "transMAX"  and  "transMIN"  regions  receive  contributions  from  the  beam-beam 
remnants  (see  Fig.  1-1).  Thus  the  "transMIN"  region  is  very  sensitive  to  the  beam-beam 
remnants,  while  the  "transMAX"  minus  the  "transMIN"  is  very  sensitive  to  initial  and 
final-state  radiation. 


6.1  The  MAX/MIN  Transverse  Regions 
As  shown  in  Figure  6-3  we  use  the  direction  of  the  highest  Pt  jet  in  the  region  |r||  < 
2,  jet#l,  to  define  the  two  "transverse"  regions,  60°  <  |A(t)|  <  120°  and  60°  <  -\A^  \  <  120°. 
On  an  event-by-event  basis,  we  define  "transMAX"  and  "transMIN"  to  be  the  maximum 
and  minimum  of  these  two  regions.  "TransMAX"  and  "transMIN"  each  have  an  area  in 
r\-^  space  of  Ar\A^  =  4nJ6.  When  looking  at  multiplicities  MAX  and  MDSf  refer  to  the 
number  of  charged  particles.  When  we  consider  Pjsum,  then  MAX  and  MIN  refer  to  the 
scalar  px  sum  of  charged  particles  and  when  we  consider  Ersum,  then  MAX  and  MIN 
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refer  to  the  scalar  Et  sum  of  particles  (or  calorimeter  towers).  The  overall  "transverse" 
region  corresponds  to  the  average  of  the  "transMAX"  and  "transMIN"  densities. 

As  illustrated  in  Fig.  6-4,  one  expects  that  "transMAX"  will  pick  up  the  hardest 
initial  or  final-state  radiation  while  both  "transMAX"  and  "transMIN"  should  receive 
"beam-beam  remnant"  contributions.  Hence  one  expects  "transMIN"  to  be  more 
sensitive  to  the  "beam-beam  remnant"  component  of  the  "underlying  event".  This  idea, 
was  first  suggested  by  Bryan  Webber,  and  implemented  by  in  a  paper  by  Jon  Pumplin 
[7,51-53].  Also,  Valaria  Tano  [54,  55]  studied  this  in  her  Run  1  analysis  of  maximum 
and  minimum  transverse  cones  (R  =  0.7). 

6.2  "Leading  Jet"  Events 

Figures  6-5  thru  6-13  show  the  data  on  the  observables  in  Table  4-3  at  1.96  TeV  for 
"leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with 
PYTHIA  Tune  A  and  HERWIG.  The  plots  show  the  uncorrected  data  (with  statistical 
errors  only)  compared  with  the  theory  after  detector  simulation  (CDFSIM).  The  plots 
also  shows  the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the 
statistical  error  and  the  systematic  uncertainty  as  described  in  Table  5-4)  compared  with 
the  theory  at  the  particle  level  (i.e.  generator  level). 

Figures  6-5  thru  6-8  reflect  a  common  trend  in  the  data.  HERWIG  is  consistently 
below  the  data  and  PYTHIA  Tune  A  for  leading  jet  Pj  less  than  about  150  GeV.  It  is 
interesting,  however,  that  HERWIG  agrees  well  for  PT(jet#l)  >  150  GeV  for  the  average 
density  of  charged  particles  and  average  charged  PTsum  density.  The  "transMIN" 
densities  are  more  sensitive  to  the  "beam-beam  remnant"  and  multiple  parton  interaction 
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components  of  the  "underlying  event"  and  PYHTIA  Tune  A  (with  multiple  parton 
interactions)  does  a  better  job  describing  the  data  than  HER  WIG  (without  multiple  parton 
interactions). 


"TransMAX/MIN"  Charged  Particle  Density:  clN/diid(|) 


1.4 


S-1.2 


CDF  Run  2  Preliminary 

data  uncorrected 

theory  +  COFSIM 


100  150  200  250  300  350 

PT(jet#1  uncorrected)  (GeV/c) 


450 


TranslUIAX/IUIIN"  Cliarged  Particle  Density:  dN/diid(|> 


1.4 
^1.2 

« 

c 

2  1.0 

I 


"transMAX" 


CDF  Run  2  Preliminary 

'  ^9ta  corrected  to  particle  levet  .^^    m^»    *^- ^ 


r  0.2 


150    200    250    300    350 
PT(jet#1)  (GeV/c) 


400 


450 


Figure  6-5:  Data  at  1.96  TeV  on  the  density  of  charged  particles,  dNchg/dridc}),  with  px  > 
0.5  GeV/c  and  |ti|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with  PYTHIA 
Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 
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Figure  6-6:  Data  at  1 .96  TeV  on  the  density  of  charged  particles,  dNchg/dridcj),  with  pr  > 
0.5  GeV/c  and  |ti|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj 
compared  with  PYTHIA  Tune  A  and  HER  WIG.  {top)  Shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation  (CDFSIM). 
(bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the 
statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the  particle 
level  (i.e.  generator  level). 
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Figure  6-7:  Data  at  1 .96  TeV  on  the  PTsum  density  of  charged  particles,  dPTsum/dT|d(t), 
with  pt  >  0.5  GeV/c  and  |t||  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "leading 
jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with  PYTHIA 
Tune  A  and  HERWIG.  {top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  {bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  {i.e.  generator 
level). 
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Figure  6-8:  Data  at  1 .96  TeV  on  the  charged  PTsum  density,  dPTsum/drid(j),  with  px  >  0.5 
GeV/c  and  |ri|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and  "transMIN") 
for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Ft  compared 
with  PYTHIA  Tune  A  and  HER  WIG.  (top)  Shows  the  uncorrected  data  (with  statistical 
errors  only)  compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom) 
Shows  the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical 
error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e. 
generator  level). 


Figure  6-9  compares  PYTHIA  Tune  A  and  HERWIG  with  the  data  on  the  average 
<Pt>  of  charged  particles  in  the  "transverse"  region  for  "leading  jet"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Pj.  The  top  shows  the  uncorrected  data  (with 
statistical  errors  only)  compared  with  the  theory  after  detector  simulation  (CDFSIM), 
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while  the  bottom  shows  the  data  corrected  to  the  particle  level  (with  errors  that  include 
both  the  statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the 
particle  level  (i.e.  generator  level).  HER  WIG,  even  after  the  data  is  corrected,  predicts  a 
softer  average  <pT>  of  charged  particles  in  the  "transverse"  region  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj. 
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Figure  6-9:  Data  at  1 .96  TeV  on  the  average  <Pt>  of  charged  particles  with  px  >  0.5 
GeV/c  and  |r||  <  1  in  the  "transverse"  region  for  "leading  jet"  events  defined  in  Fig.  6-3  as 
a  function  of  the  leading  jet  Pj  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top) 
Shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 
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Figure  6-10  compares  PYTHIA  Tune  A  and  HERWIG  with  the  data  on  the  average 
maximum  pj,  PTmax,  for  charged  particles  in  the  "transverse"  region  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj.  The  top  shows  the 
uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after  detector 
simulation  (CDFSIM),  while  the  bottom  shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level).  The  agreement  is 
better  than  for  the  average  <Pt>,  but  HERWIG  is  still  too  soft  for  leading  jet  Pj  less  than 
about  150  GeV. 

Figures  6-11  and  6-12  compare  PYTHIA  Tune  A  and  HERWIG  with  the  data  on 
the  the  ETsum  density,  dEx/dTidc]),  for  particles  with  |ri|  <  1  in  the  transverse  region  for 
"leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px.  The  top  shows 
the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after  detector 
simulation  (CDFSIM),  while  the  bottom  shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level).  Both  PYTHIA  Tune 
A  and  HERWIG  fall  below  the  data  throughout  the  energy  range. 

Figure  6-13  compares  PYTHIA  Tune  A  and  HERWIG  with  the  data  at  on  the 
charged  fraction,  PTsum/ETsum,  in  the  "transverse"  region  for  "leading  jet"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj,  where  PTsum  includes  charged 
particles  with  px  >  0.5  GeV/c  and  |ri|  <  1  and  the  ETsum  includes  all  particles  with  |r||  < 
1.  The  top  shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the 
theory  after  detector  simulation  (CDFSIM),  while  the  bottom  shows  the  data  corrected  to 
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the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the  systematic 
uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator  level).  Both 
MC  generators  predict  a  larger  charged  fraction  in  the  transverse  region  than  is  indicated 
by  the  data. 
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Figure  6-10:  Data  at  1.96  TeV  on  the  average  maximum  pj,  PTmax,  for  charged  particles 
with  pt  >  0.5  GeV/c  and  |r||  <  1  in  the  "transverse"  region  for  "leading  jet"  events  defined 
in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with  PYTHIA  Tune  A  and 
HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with 
the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the 
particle  level  (with  errors  that  include  both  the  statistical  error  and  the  systematic 
uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 
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Figure  6-11:  Data  at  1.96  TeV  on  the  ETsum  density,  dET/dT|d(t),  for  particles  with  |t||  <  1 
in  the  "transMAX"  and  "transMIN"  regions  for  "leading  jet"  events  defined  in  Fig.  6-3  as 
a  function  of  the  leading  jet  Pj  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top) 
Shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 
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Figure  6-12:  Data  at  1.96  TeV  on  the  ETsum  density,  dEx/dTidct),  for  particles  with  |ri|  <  1 
in  the  "transverse"  region  (average  of  "transMAX"  and  "transMESf")  for  "leading  jet" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with  PYTHIA 
Tune  A  and  HER  WIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 
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Figure  6-13:  Data  at  1.96  TeV  on  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  for  "leading  jet"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Pj,  where  PTsum  includes  charged  particles  with  px  >  0.5  GeV/c  and  hi  <  1 
and  the  ETsum  includes  all  particles  with  |ri|  <  1.  The  data  are  compared  with  PYTHIA 
Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  {i.e.  generator 
level). 


6.3  "Back-to-Back"  Events 
Figure  6-14  thru  6-22  show  the  data  on  the  observables  in  Table  4-3  at  1.96  TeV  for 
"back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared 
with  PYTHIA  Tune  A  and  HERWIG.     The  plots  show  the  uncorrected  data  (with 
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statistical  errors  only)  compared  with  the  theory  after  detector  simulation  (CDFSIM). 
The  plots  also  shows  the  data  corrected  to  the  particle  level  (with  errors  that  include  both 
the  statistical  error  and  the  systematic  uncertainty  as  described  in  Table  5-4)  compared 
with  the  theory  at  the  particle  level  {i.e.  generator  level). 

Figures  6-14  thru  6-17  continue  to  exhibit  the  trend  that  HER  WIG  is  consistently 
below  the  data  and  PYTHIA  Tune  A  for  leading  jet  Pt  less  than  about  150  GeV. 
However,  the  back-to-back  data  reveal  that  the  densities  fall  with  increasing  PxCJetttl). 
Since  the  hard  initial  and  final-state  radiation  is  suppressed  in  back-to-back  events  the 
differences  between  PYHTIA  Tune  A  and  HERWIG  in  the  "beam-beam  remnant"  and 
multiple  parton  scattering  component  of  the  "underlying  event"  are  more  evident. 
PYHTIA  Tune  A  (with  multiple  parton  interactions)  does  a  better  job  describing  the  scale 
and  slope  of  the  data  than  HERWIG  (without  multiple  parton  interactions). 

Figure  6-18  compares  PYTHIA  Tune  A  and  HERWIG  with  the  data  on  the  average 
<Pt>  of  charged  particles  with  px  >  0.5  GeV/c  and  |t||  <  1  in  the  "transverse"  region  for 
"back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj.  The  top 
shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM),  while  the  bottom  shows  the  data  corrected  to  the  particle 
level  (with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  {i.e.  generator  level).  Both  PYTHIA  Tune 
A  and  HERWIG  predict  a  softer  average  <pT>  of  charged  particles  than  is  exhibited  by 
the  data. 
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Figure  6-14:  Data  at  1.96  TeV  on  the  density  of  charged  particles,  dNchg/dridcj),  with  px 
>  0.5  GeV/c  and  |ti|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "back-to-back" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px  compared  with  PYTHIA 
Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 
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Figure  6-15:  Data  at  1.96  TeV  on  the  density  of  charged  particles,  dNchg/dTid(j),  with  pj 
>  0.5  GeV/c  and  |t||  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet 
Pt  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data 
(with  statistical  errors  only)  compared  with  the  theory  after  detector  simulation 
(CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level  (with  errors  that 
include  both  the  statistical  error  and  the  systematic  uncertainty)  compared  with  the  theory 
at  the  particle  level  (i.e.  generator  level). 
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Figure  6-16:  Data  at  1.96  TeV  on  the  PTsum  density  of  charged  particles,  dPTsum/dridcl), 
with  pT  >  0.5  GeV/c  and  |ri|  <  1  in  the  "transMAX"  and  "transMIN"  regions  for  "back-to- 
back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with 
PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors 
only)  compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the 
data  corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and 
the  systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 
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Figure  6-17:  Data  at  1.96  TeV  on  the  PTsum  density,  dPTsum/dridct),  with  pj  >  0.5  GeV/c 
and  |ri|  <  1  in  the  "transverse"  region  (average  of  "transMAX"  and  "transMIN")  for 
"back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Ft  compared 
with  PYTHIA  Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical 
errors  only)  compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom) 
Shows  the  data  corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical 
error  and  the  systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e. 
generator  level). 


93 


"Transverse"  Charged  Particle  Mean  PT 


2.0 

> 
t 

A 

1 


2  10 

I 


0.5 


CDF  Run  2  Preliminary 

data  uncorrected 
thaory  -f  CDFSIM 

I^BTeV 


"Back-to-Back" 


PY  Tune  A 


^ i 


-^-'f-'- 


Midpoint  R  =  0.7  |iiaet#1)  <  2 
Charged  Particles  (|7i|<1.0,  PT>0.5  GeV/c) 


100         150         200         250         300         350 
PT(jet#1  uncorrected)  (GeV/c) 


400 


450 


"Transverse"  Charged  Particle  Mean  PT 


2.0 


> 
o 
(3 


0-   1.5 


1 
a 


8  1-0 
I 

c 
0.5 


CDF  Run  2  Preliminary 

data  corrected  to  particle  level 

1.96  TeV 


"Back-to-Back" 


5} 


gmmiiiitiiHiHi  i  iVl  i  i  J  ^j 

*  V  J- 


HW 


IMidPoint  R  =  0.7  |Tl(jet#1)  <  2 

Charged  Particles  (|t||<1.0,  PT>0.5  GeV/c) 

1 1 \ 1 


50     100    150    200    250    300 
PT(jet#1)  (GeV/c) 


350    400 


450 


Figure  6-18:  Data  at  1.96  TeV  on  the  average  <pt>  of  charged  particles  with  px  >  0.5 
GeV/c  and  |ti|  <  1  in  the  "transverse"  region  for  "back-to-back"  events  defined  in  Fig.  6-3 
as  a  function  of  the  leading  jet  Px  compared  with  PYTHIA  Tune  A  and  HERWIG.  (top) 
Shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 

Figure  6-19  compares  PYTHIA  Tune  A  and  HERWIG  with  the  data  on  the  average 

maximum  px,  PTmax,  for  charged  particles  with  px  >  0.5  GeV/c  and  |ti|  <  1  in  the 

"transverse"  region  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the 

leading  jet  Pt-  The  top  shows  the  uncorrected  data  (with  statistical  errors  only)  compared 

with  the  theory  after  detector  simulation  (CDFSIM),  while  the  bottom  shows  the  data 
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corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level).  Both  MC  generators  predict  a  too  soft  average  maximum  pr,  but  only  PYTHIA 
Tune  A  correctly  reproduces  the  decreasing  slope  of  the  observable. 
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Figure  6-19:  Data  at  1.96  TeV  on  the  average  maximum  px,  PTmax,  for  charged  particles 
with  Px  >  0.5  GeV/c  and  |r||  <  1  in  the  "transverse"  region  for  "back-to-back"  events 
defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px  compared  with  PYTHIA  Tune  A 
and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only)  compared 
with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected 
to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the  systematic 
uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 
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Figures  6-20  and  6-21  compare  PYTHIA  Tune  A  and  HER  WIG  with  the  data  on 
the  ETsum  density,  dEx/dridcl),  for  particles  with  |ti|  <  1  in  the  transverse  region  for 
"back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt-  The  top 
shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM),  while  the  bottom  shows  the  data  corrected  to  the  particle 
level  (with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level).  Both  MC  generators 
predict  too  little  ETsum  density,  but  only  PYTHIA  Tune  A  correctly  reproduces  the 
decreasing  slope  of  the  observable. 

Figure  6-22  compares  PYTHIA  Tune  A  and  HER  WIG  with  the  data  on  the  charged 
fraction,  PTsum/ETsum,  in  the  "transverse"  region  for  "back-to-back"  events  defined  in 
Fig.  6-3  as  a  function  of  the  leading  jet  Pt,  where  PTsum  includes  charged  particles  with 
Pt  >  0.5  GeV/c  and  h|  <  1  and  the  ETsum  includes  all  particles  with  |ti|  <  1.  The  top 
shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM),  while  the  bottom  shows  the  data  corrected  to  the  particle 
level  (with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level).  Both  MC  generators 
predict  a  larger  charged  fraction  in  the  transverse  region  than  is  indicated  by  the  data. 
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Figure  6-20:  Data  at  1.96  TeV  on  the  ETsum  density,  dEx/dridct),  for  particles  with  hi  <  1 
in  the  "transMAX"  and  "transMIN"  regions  for  "back-to-back"  events  defined  in  Fig.  6-3 
as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA  Tune  A  and  HER  WIG.  (top) 
Shows  the  uncorrected  data  (with  statistical  errors  only)  compared  with  the  theory  after 
detector  simulation  (CDFSIM).  (bottom)  Shows  the  data  corrected  to  the  particle  level 
(with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty) 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 
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Figure  6-21:  Data  at  1.96  TeV  on  the  ETsum  density,  dEj/dTidct),  for  particles  with  |ri|  <  1 
in  the  "transverse"  region  (average  of  "transMAX"  and  "transMIN")  for  "back-to-back" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA 
Tune  A  and  HERWIG.  (top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 
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Figure  6-22:  Data  at  1.96  TeV  on  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  for  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the 
leading  jet  Pj,  where  PTsum  includes  charged  particles  with  px  >  0.5  GeV/c  and  hi  <  1 
and  the  ETsum  includes  all  particles  with  |ri|  <  1.  The  data  are  compared  with  PYTHIA 
Tune  A  and  HERWIG.  {top)  Shows  the  uncorrected  data  (with  statistical  errors  only) 
compared  with  the  theory  after  detector  simulation  (CDFSIM).  (bottom)  Shows  the  data 
corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 


6.4  "Leading  Jet"  versus  "Back-to-Back"  Events 
Figure  6-23  and  Figure  6-24  compare  the  data  on  the  density  of  charged  particles 
and  the  charged  PTsum  density  in  the  "transverse"  region  corrected  to  the  particle  level 
for  "leading  jet"  and  "back-to-back"  events  with  PYTHIA  Tune  A  and  HERWIG  at  the 
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particle  level.  As  expected,  the  "leading  jet"  and  "back-to-back"  events  behave  quite 
differently.  For  the  "leading  jet"  case  the  "transMAX"  densities  rise  with  increasing 
PxCJetttl),  while  for  the  "back-to-back"  case  they  fall  with  increasing  PxGetttl).  The  rise 
in  the  "leading  jet"  case  is,  of  course,  due  to  hard  initial  and  final-state  radiation,  which 
has  been  suppressed  in  the  "back-to-back"  events.  The  "back-to-back"  events  allow  for 
closer  look  at  the  "beam-beam  remnant"  and  multiple  parton  scattering  component  of  the 
"underlying  event"  and  PYTHIA  Tune  A  (with  multiple  parton  interactions)  does  a  better 
job  describing  the  data  than  HER  WIG  (without  multiple  parton  interactions). 

The  "transMIN"  densities  are  more  sensitive  to  the  "beam-beam  remnant"  and 
multiple  parton  interaction  component  of  the  "underlying  event".  The  "back-to-back" 
data  show  a  decrease  in  the  "transME^"  densities  with  increasing  PrOetttl)  which  is 
described  fairly  well  by  PYTHIA  Tune  A  (with  multiple  parton  interactions)  but  not  by 
HERWIG  (without  multiple  parton  interactions).  The  decrease  of  the  "transMIN" 
densities  with  increasing  PjCJetttl)  for  the  "back-to-back"  events  is  very  interesting  and 
might  be  due  to  a  "saturation"  of  the  multiple  parton  interactions  at  small  impact 
parameter.  Such  an  effect  is  included  in  PYTHIA  Tune  A  but  not  in  HERWIG  (without 
multiple  parton  interactions). 


100 


"TransMAX"  Charged  Particle  Density:  dN/dTid(t» 


50     100    150    200    250    300    350    400    450 
PT(jet#1)  (GeV/c) 


"TranslUIIN"  Charged  Particle  Density:  dN/dridtj) 


0.6 


0.5 


CDF  Run  2  Preliminary 

data  correctedto  particle  levei- 
1.96  TeV 


MIdPoint  R  =  0.7  haet#1)  <  2 


0  50  100  150  200  250  300  350  400  450 

PT(jet#1)  (GeV/c) 


"Transverse"  Charged  Particle  Density:  dN/dT|di|) 


1.0 


CDF  Run  2  Preliminary 

data  corrected  to  particle  level 


Charged  Particles  (hll<10,  PT>0.5  GeV/c) 

H 1 1 \ \- 

200    250    300    350    400    450 
PT(jet#1)  (GeV/c) 


Figure  6-23:  Data  at  1 .96  TeV  on  the  density  of  charged  particles,  dNchg/dr|d(t),  with  px  > 
0.5  GeV/c  and  |ri|  <  1  in  the  "transMAX"  region  (top),  "transMIN"  region  (middle),  and 
"transverse"  region  (average  of  "transMAX"  and  "transMIN")  (bottom)  for  "leading  jet" 
and  "back-to-back"  events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Px 
compared  with  PYTHIA  Tune  A  and  HER  WIG.  The  data  are  corrected  to  the  particle 
level  (with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty)  and 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 
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Figure  6-24:  Data  at  1.96  TeV  on  charged  PTsum  density  of  charged  particles, 
dPTsum/dr|d(l),  with  px  >  0.5  GeV/c  and  |ti|  <  1  in  the  "transMAX"  region  (top), 
"transMIN"  region  (middle),  and  "transverse"  region  (average  of  "transMAX"  and 
"transMIN")  (bottom)  for  "leading  jet"  and  "back-to-back"  events  defined  in  Fig.  6-3  as  a 
function  of  the  leading  jet  Px  compared  with  PYTHIA  Tune  A  and  HERWIG.  The  data 
are  corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical  error  and  the 
systematic  uncertainty)  and  compared  with  the  theory  at  the  particle  level  (i.e.  generator 
level). 
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Fig.  6-25  compares  the  data  on  average  <Pt>  of  charged  particles  and  the  average 
maximum  charge  particle  px,  PTmax,  in  the  "transverse"  region  corrected  to  the  particle 
level  for  "leading  jet"  and  "back-to-back"  events  with  PYTHIA  Tune  A  and  HER  WIG  at 
the  particle  level.  Again  the  "leading  jet"  and  "back-to-back"  events  behave  quite 
differently. 

Fig.  6-26  shows  the  data  corrected  to  the  particle  level  on  the  ETsum  density, 
dEi/dridct),  in  the  "transverse"  region  for  "leading  jet"  and  "back-to-back"  events 
compared  with  PYTHIA  Tune  A  and  HERWIG  at  the  particle  level.  Neither  PYTHIA 
Tune  A  or  HERWIG  produce  enough  energy  in  the  "transverse"  region.  HERWIG  has 
more  "soft"  particles  than  PYTHIA  Tune  A  does  slightly  better  in  describing  the  energy 
density  in  the  "transMAX"  and  "transMIN"  region. 

Fig.  6-27  shows  the  data  corrected  on  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  for  "leading  jet"  and  "back-to-back"  events  compared  with  PYTHIA 
Tune  A  and  HERWIG  at  the  particle  level.  Neither  PYTHIA  Tune  A  or  HERWIG 
produce  enough  energy  in  the  "transverse"  region  and  therefore  predict  too  large  of  a 
charged  fraction.  Note  that  both  PYTHIA  Tune  A  and  HERWIG  predict  a  charged 
fraction  of  about  0.5  if  one  includes  all  particles  in  both  Pjsum  and  Ejsum. 
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Figure  6-25:  Data  at  1.96  TeV  on  the  average  <pT>  of  charged  particles  (top)  and  the 
average  maximum  pj,  PTmax,  for  charged  particles  (bottom)  with  pj  >  0.5  GeV/c  and  |ri| 
<  1  in  the  "transverse"  region  for  "leading  jet"  and  "back-to-back"  events  defined  in  Fig. 
6-3  as  a  function  of  the  leading  jet  Pt  compared  with  PYTHIA  Tune  A  and  HER  WIG. 
The  data  are  corrected  to  the  particle  level  (with  errors  that  include  both  the  statistical 
error  and  the  systematic  uncertainty)  and  compared  with  the  theory  at  the  particle  level 
(i.e.  generator  level). 
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Figure  6-26:  Data  at  1.96  TeV  on  the  ETsum  density,  dET/dT|d(|),  for  particles  with  |ti|  <  1 
in  the  "transMAX"  region  (top),  "transMEN"  region  (middle),  and  "transverse"  region 
(average  of  "transMAX"  and  "transMEN")  (bottom)  for  "leading  jet"  and  "back-to-back" 
events  defined  in  Fig.  6-3  as  a  function  of  the  leading  jet  Pj  compared  with  PYTHIA 
Tune  A  and  HERWIG.  The  data  are  corrected  to  the  particle  level  (with  errors  that 
include  both  the  statistical  error  and  the  systematic  uncertainty)  and  compared  with  the 
theory  at  the  particle  level  (i.e.  generator  level). 
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Figure  6-27:  Data  at  1.96  TeV  on  the  charged  fraction,  PTsum/ETsum,  in  the 
"transverse"  region  defined  in  Fig.  6-3  for  "leading  jet"  events  (top)  and  "back-to-back" 
events  (bottom)  as  a  function  of  the  leading  jet  Pj,  where  PTsum  includes  charged 
particles  with  px  >  0.5  GeV/c  and  hi  <  1  and  the  ETsum  includes  all  particles  with  |ti|  < 
1,  compared  with  PYTHIA  Tune  A  and  HERWIG.  The  data  are  corrected  to  the  particle 
level  (with  errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty)  and 
compared  with  the  theory  at  the  particle  level  (i.e.  generator  level). 

Figure  6-28  compares  the  corrected  "leading  jet"  and  "back-to-back"  data  on  the 
density  of  charged  particles,  dN/dT|d(t),  the  Pisum  density,  dPrsum/  dridc]),  and  the  Exsum 
density,  dExsun/dridcl),  for  "TransDIF'  (i.e.  the  "transMAX"  region  minus  the  "transMIN" 
region)  as  a  function  of  the  leading  jet  Px.  Figure  6-28  also  shows  the  predictions  of 
PYTHIA  Tune  A  and  HERWIG.   The  difference  between  "transMAX"  and  "transMIN" 
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is  very  sensitive  to  the  "hard  scattering"  component  of  the  "underiying  event"  (i.e  hard 
initial  and  final-state  radiation)  and  therefore  we  expect  to  see  a  big  difference  between 
"leading  jet"  events  and  the  "back-to-back"  events  which  is  the  case.  Again,  PYTHIA 
Tune  A  does  a  better  job  describing  the  data  than  does  HERWIG  fro  the  charged  particle 
densities  and  the  charged  Ptsuhi  densities.  Surprisingly,  both  PYTHIA  Tune  A  and 
HERWIG  agree  with  the  data  for  the  "transDIF'  Exsum  density. 
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Figure  6-28:  Data  at  1.96  TeV  on  the  density  of  charged  particles,  dN/dridc])  and  the  Pxsum 
density,  dPrsum/  dridcj)  for  charged  particles  with  pr  >  0.5  GeV/c  and  hi  <  1  for  the 
"transMAX"  minus  the  "transMIN"  region  defined  in  figure  6-3  for  "leading  jet"  and 
"back-to-back"  events  defined  in  Fig.  6-2  as  a  function  of  the  leading  jet  Pt  compared 
with  PYTHIA  Tune  A  and  HERWIG.  The  "transDIF'  Exsum  density,  dExsum/dTidct)  for 
particles  with  |ti|  <  1  is  also  presented.  The  data  are  corrected  to  the  particle  level  (with 
errors  that  include  both  the  statistical  error  and  the  systematic  uncertainty)  and  compared 
with  the  theory  at  the  particle  level  (i.e.  generator  level). 


CHAPTER  7 
SUMMARY  AND  CONCLUSION 

The  goal  of  this  analysis  is  to  produce  data  on  the  "underlying  event"  that  is 
corrected  to  the  particle  level  so  that  it  can  be  used  to  tune  the  QCD  Monte-Carlo  models 
without  requiring  CDF  detector  simulation  {i.e.  CDFSIM).  Unlike  our  previous  Run  2 
"underlying  event"  analysis  [6]  which  used  JetClu  to  define  "jets"  and  compared 
uncorrected  data  with  PYTHIA  Tune  A  and  HERWIG  after  detector  simulation 
(CDFSIM),  in  this  analysis  we  use  the  MidPoint  algorithm  (R  =  0.7,  fmerge  =  0.75)  and 
correct  the  observables  to  the  particle  level.  The  corrected  observables  are  then 
compared  with  PYTHIA  Tune  A  and  HERWIG  at  the  particle  level  (i.e.  generator  level). 
This  analysis  also  provides  our  first  look  at  the  energy  in  the  "underlying  event"  (i.e.  the 
"transverse"  region). 

In  this  analysis  we  look  at  both  the  charged  particle  and  the  energy  components  of 
the  "underlying  event".  We  use  the  direction  of  the  leading  calorimeter  jet  in  each  event 
to  define  two  "transverse"  regions  of  r\-(^  space  that  are  very  sensitive  to  the  "underlying 
event".  In  addition,  by  selecting  events  with  at  least  two  jets  that  are  nearly  back-to-back 
(A(t)i2  >  150°)  with  PT(jet#3)  <  15  GeV/c  we  are  able  to  look  closer  at  the  "beam-beam 
remnant"  and  multiple  parton  interaction  components  of  the  "underlying  event". 

Comparing  the  corrected  observables  with  PYTHIA  Tune  A  and  HERWIG  at  the 
particle  level  (i.e.  generator  level)  continue  to  validate  the  conclusions  we  found  when 
comparing  the  uncorrected  data  with  the  Monte-Carlo  models  after  detector  simulation 
(i.e.  CDFSIM)  [8].  PYTHIA  Tune  A  (with  multiple  parton  interactions)  does  a  better  job 
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in  describing  the  "underiying  event"  (i.e.  "transverse"  regions)  for  both  "leading  jet"  and 
"back-to-back"  events  than  does  HER  WIG  (without  multiple  parton  interactions). 
Herwig  does  not  have  enough  activity  in  the  "underlying  event"  for  PrCJetttl)  less  than 
about  150  GeV,  which  was  also  observed  in  a  published  Run  I  analysis  [13]. 

Figure  7-1  compares  the  original  Run  I  published  data  on  the  underlying  event  with 
our  current  results.  We  have  extended  the  range  of  the  leading  jet  Pj  from  50  GeV  to 
over  400  GeV.  We  now  use  the  MidPoint  Jet  finding  algorithm  and  most  importantly 
correct  the  data  back  to  the  particle  level. 

This  analysis  gives  our  first  look  at  the  energy  in  the  "underlying  event"  (i.e.  the 
"transverse"  region).  Figure  6-28  shows  that  PYTHIA  Tune  A  and  HERWIG  agree  with 
the  "transDIF'  Eisum  density  data.  Since  both  MC  generators  correctly  predict  the  "hard 
scattering"  component  of  the  underlying  event,  this  strongly  suggests  that  the 
discrepancies  and  low  jet#l  Pt  reflect  a  needed  tuning  of  the  "beam-beam  remnant"  and 
multiple  parton  interaction  components  of  the  underlying  event. 
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Figure  7-1:  A  comparison  between  the  Run  I  and  Run  11  analysis  on  the  PTsum  density  in 
the  "transverse"  region,  (top)  Data  uncorrected  (with  errors  that  include  both  statistical 
and  the  systematic  uncertainty)  at  1.8  TeV  on  the  charged  PTsum  density,  dPTsum/dr|d(|), 
with  PT  >  0.5  GeV/c  and  |t||  <  1.0  for  leading  jet  events  defined  in  Fig.  6-3  as  a  function 
of  leading  jet  PT  compared  PYTHIA  set  A,  PYTHIA  set  B,  and  HERWIG.  (bottom)  Data 
corrected  (with  errors  that  include  both  statistical  and  the  systematic  uncertainty)  at  1 .96 
TeV  on  the  charged  PTsum  density,  dPTsura/dTid(t),  with  PT  >  0.5  GeV/c  and  |ri|  <  1.0  for 
"leading  jet  events"  defined  in  Fig.  6-3  as  a  function  of  leading  jet  PT  compared  PYTHIA 
Tune  A  and  HERWIG. 
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Elementary  particle  physics  attempts  to  answer  very  fundamental  questions  of  how 

the  Universe  was  created  and  how  it  works.   One  approach  to  unlocking  these  mysteries 

is  by  colliding  very  fast  moving  protons  and  antiprotons  and  studying  the  outcome.  Most 

of  the  time  these  particles  ooze  through  each  other,  but  occasionally  we  get  a  collision 

that  is  characterized  by  a  large  amount  of  transverse  momentum.   This  signals  a  special 

kind  of  collision  that  can  be  calculated  by  the  theorist.  The  problem  lies  in  the  fact  that 

only  a  portion  of  the  collision  can  be  calculated.   The  majority  of  the  collision  is  messy 

and  must  be  modeled.  The  data  presented  here  helps  to  improve  the  current  models  and 

allows  for  a  better  understanding  of  the  dynamics  of  nuclear  forces. 
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